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In  late  1976  a  numerical  wave  hindcast  model  was  developed  by  the  Wave  In- 
formation Study  (WIS)  at  the  U.S.  Army  Engineer  Waterways  Experiment  Station 
(WES)  and  it  was  employed  to  generate  20  years  of  wave  conditions  in  10-meter 
water  depth  along  the  eastern  seaboard  of  the  United  States.  The  model  comprises 
of  deepwater  and  shallow  water  submodels.  The  significant  nonlinear  wave- wave 
interaction  process  is  included  in  the  deepwater  submodel  based  on  the  study  of 
parametric  model.  Recently,  more  detailed  field  wave  information  became  available 
in  both  offshore  and  nearshore  regions,  particularly  shallow  water  directional  wave 
information  reduced  from  data  collected  along  the  Florida  Coast.  Based  upon  those 
data,  the  WIS  model  is  reexamined. 

A  new  wave  model,  named  here  as  the  Unified  Coastal  Wave  Prediction  (UCWP) 
model,  is  developed.  This  model  is  of  a  discrete  type  and  it  links  a  wind  field 
model  to  a  combined  deepwater  and  shallow  water  wave  prediction  model.  The 
deepwater  wave  submodel  follows  the  WIS  model  but  with  revised  source  terms 
and  dissipation  terms.  The  shallow  water  wave  submodel  utilizies  the  refraction 
law  and  the  equation  of  conservation  of  energy  to  solve  for  the  shallow  water  waves. 
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The  UCWP  model  represents  a  significant  improvement  over  existing  models, 
at  least  for  applications  to  the  southeastern  region  of  the  United  States.  The  im- 
provement is  particularly  evident  in  wave  period  and  wave  direction.  In  addition  to 
the  development  of  the  new  numerical  model,  a  number  of  fundamental  aspects  of 
air-sea  interactions  are  also  examined  in  light  of  the  field  data.  New  semi-empirical 
formulas  are  proposed  for  a  number  of  importajit  parameters,  including  the  JON- 
SWAP  spectral  parameters,  required  to  define  wind  wave  spectra. 


CHAPTER  1 
INTRODUCTION 


1.1     Background 


Wave  information  is  required  for  a  variety  of  purposes,  including  design  of 
coastal  ajid  offshore  structures,  planning  of  marine  operations,  ship  design,  coastal 
process  studies,  and  coastal  erosion  prevention,  etc.  Because  of  a  general  lack  of 
measured  wave  data,  particularly  long-term  data,  wave  hindcast,  which  is  a  tech- 
nique in  calculating  the  events  of  wind-generated  water  surface  waves  in  the  past,  is 
often  the  best  means  to  estimate  representative  wave  statistics  and  is  almost  always 
the  only  meajis  to  estimate  characteristics  of  the  infrequently  occurring  extreme 
waves. 

In  early  days,  methods  utilized  for  hindcasting  wind  waves  involved  empirically 
based  nomographys  and  tables  to  evaluate  wave  height  and  wave  period  for  promi- 
nent waves.  A  well-known  example  is  the  Sverdrup-Munk-Bretschneider  (SMB) 
model  which  was  established  by  matching  some  predescribed  math  functions  to 
measured  wave  heights  and  wave  periods.  The  model  was  established  in  a  period 
during  and  after  World  War  II  and  was  used  exclusively  for  either  hindcast  or  fore- 
cast of  ocean  waves  before  more  modern  wave  models  were  developed  later  in  the 
1970s.  With  more  physical  processes  governing  the  nature  of  waves  being  known 
and  concepts  of  wave  energy  spectra  becoming  widely  accepted,  spectral  models 
were  then  developed  for  wave  hindcast,  or  forecast. 

The  spectral  models  are  commonly  categorized  into  two  major  branches:  the 
discrete  model  and  the  parametric  model.  With  wind  information  as  input,  both 
the  discrete  model  and  the  parametric  model  generate  wave  energy  as  a  function 
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of  frequency  and  direction,  propagate  energy  across  oceanic  region,  and  consider 
energy  spreading  ajid  dissipation.  The  difference  between  the  two  models  is  that 
the  discrete  model  evaluates  wave  spectra  by  summing  all  individual  spectral  com- 
ponents corresponding  to  different  physical  processes,  while  the  parametric  model 
constructs  wave  spectra  by  solving  numerically  spectral  parameters  from  the  prog- 
nostic equations  which  represent  the  projections  of  wind-wave  relationships  in  the 
spectral  parameter  space. 

The  advantage  of  using  parametric  model  to  hindcast  the  wind  wave  spectra  is 
that  its  numerical  solution  procedure  is  simple.  Computations  in  solving  spectral 
pcirELmeters  in  the  parametric  model  axe  solely  performed  in  the  frequency  domain. 
The  directional  spreading  of  wave  energy  is  considered  in  the  model  by  multiplying 
the  resulting  frequency  spectrum  with  an  angular  distribution  function.  The  as- 
sumed angular  distribution  function  is  a  normalized  bell-shaped,  symmetrical  func- 
tion with  respect  to  local  wind  direction.  The  disadvantage  of  using  the  parametric 
model  is  that  it  does  not  include  swell  and,  hence,  cannot  rigorously  consider  all  the 
physical  processes  involved  in  wave  transformation.  Furthermore,  the  parametric 
model  can  not  properly  generate  wave  spectra  in  shallow  water  where  bottom  effects 
become  important.  The  discrete  model  in  more  flexible  and  can  include  swells  just 
as  easily  as  that  of  the  regular  wind  waves.  It  can  also  handle  easily  the  shallow 
water  effects,  such  as  refrziction,  shoaling,  bottom  friction,  etc.  However,  in  the 
discrete  model  computations  involved  are  usually  tedious  and  time-consuming  since 
they  have  to  be  carried  out  in  the  combined  frequency  and  direction  domain. 

Currently,  there  are  a  number  of  spectral  models,  both  discrete  aind  parametric 
types,  available  for  wave  hindcasting.  Most  of  them  are  for  deepwater  wave  applica- 
tion. These  wave  models  are  based  upon  similar  theoretical  considerations  but  their 
numerical  methods  and  computational  requirements  are  vastly  different.  They  all 
contain  empirical  functions  or  coefficients  derived  from  actual  wave  measurements 


3 

and  aje  calibrated  on  limited  data  base.  Due  to  the  scarcity  of  the  measured  wave 
data,  none  of  the  present  models  were  tested  against  a  wide  rajige  of  wind  and  wave 
conditions.  Compaxisons  for  extreme  waves  during  severe  storms  and  hiirricanes 
were  few  and  compairisons  with  measured  wave  energy  directional  spectra  were  ex- 
tremely rare.  Accordingly,  there  is  considerable  uncertainty  about  the  accxiracy  of 
the  existing  wave  models. 

A  few  comparisons  of  different  wave  m.odels  were  made  with  each  other  in  the 
past.  Among  them,  the  work  by  Resio  and  Vincent  (1979)  offered  an  interesting 
comparison  among  the  SMB  model,  Hasselmann's  (1976)  parametric  model,  a  ver- 
sion of  the  Pierson  et  al.(1966)  model,  Barnett's  (1968)  model,  and  a  discrete  model 
by  Resio  and  Vincent  (1977,  1978).  Comparisons  of  wave  height  versus  fetch  and 
duration  for  selected  constant  wind  speeds  indicate  significant  differences  among 
these  models  and,  thus,  raise  concern  for  their  accuracy. 

There  axe  several  problems  which  need  to  be  addressed  if  wave  models  are  to 
improve  their  accuracy  and  to  be  applicable  for  a  wide  range  of  conditions  in  both 
deep  ajid  shallow  waters.  First,  source  and  sink  terms  should  be  continuously  re- 
fined with  availability  of  more  and  better  wave  data.  Second,  wave  models  should 
be  calibrated  with  directional  wave  information  and  with  a  wide  range  of  wind  con- 
ditions, i.e.,  large  systems,  small  systems,  extreme  events,  etc.  Finally,  there  is  a 
need  to  acquire  the  accurate  wind  information  since  wave  model  relies  on  wind  field 
as  input.  The  wind  information  can  be  either  measured  from  the  field  or  generated 
by  a  wind  model,  based  on  the  sea  surface  pressure  distribution.  Accurate  specifi- 
cation of  wind  fields  through  a  numerical  model  is  sometimes  difficult.  Therefore,  it 
is  necessary  to  improve  the  wind  model  or  to  properly  incorporate  measured  wind 
data  as  input  to  the  wave  model. 
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1.2     Scope  of  Work 

Recently,  more  detailed  field  wave  information  became  available  in  both  deepwa- 
ter  and  shallow  water  regions.  For  instajice,  long  term  wave  data  have  been  collected 
in  deep  water  by  several  NOAA  (National  Oceamic  and  Atmospheric  Adminstration) 
maintained  buoys  in  the  Northwest  Atlantic  Ocean  and  in  shallow  water  by  twelve 
CDN  (Coastal  Data  Network)  underwater  stations  along  the  Florida  coast.  Based 
upon  these  data,  a  new  coupled  wave  hindcast  model  is  developed.  This  model, 
referred  to  here  as  the  Unified  Coastal  Wave  Prediction  (UCWP)  model,  is  of  a  dis- 
crete type.  It  links  a  wind  field  model  to  a  combined  deepwater  and  shallow  water 
wave  prediction  model.  The  deepwater  submodel  follows  the  same  technique  used 
in  WIS  with  revised  source  and  dissipation  terms.  The  shallow  water  submodel 
utilizes  the  algorithm  developed  by  Chen  and  Wang  (1983)  with  modified  source 
term  proposed  by  Vincent  (1984).  The  model  is  unified  in  the  sense  that  the  deep- 
water  model  and  the  shallow  water  model  are  based  on  the  same  set  of  fundamental 
equations;  both  employ  discrete  finite  different  scheme  with  compatible  source  and 
dissipation  terms;  the  models  axe  on  the  same  level  in  detail  and  are  expected  to 
yield  same  degree  of  accuracy. 

A  reexamination  of  two  previously  developed  deepwater  wave  models,  the  WIS 
model  and  the  parametric  model  by  Hasselmann  et  al.(1973),  has  revealed  some 
differences  about  the  source  functions  utilized  by  the  two  models.  First,  the  non- 
conservative  term  of  wave  energy  dissipation  is  completely  neglected  from  the  para- 
metric model  while  it  is  included  in  the  discrete  model  to  account  for  energy  loss 
of  low  frequency  components.  Second,  the  conservative  wave  energy  transfer,  as 
a  result  of  nonlinear  wave-wave  interactions,  is  included  in  the  parajnetric  model 
based  on  JONSWAP  spectrum  while  it  is  included  in  the  discrete  model  based  on 
Kitaigorodskii  spectrum.  Third,  the  wave  growth  due  to  wind  input  is  considered  in- 
dependent to  the  nonlinear  wave-wave  interactions  in  the  parametric  model,  while 
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it  is  assumed  to  depend  on  the  nonlinear  wave-wave  interactions  in  the  discrete 
model.  A  well-known  problem  about  the  parametric  model  is  its  inability  of  treat- 
ing swell  components.  The  parametric  is  preferred  by  some  due  to  its  simplicity  in 
computations,  while  the  discrete  model  is  preferred  by  others  because  of  its  straight 
forward  physical  interpretation. 

Inspection  of  field  wave  data  collected  by  NOAA  buoys  reveals  that  nonconser- 
vative  turbulent  dissipation  of  wave  energy  induced  by  wave  motion  should  not  be 
neglected  in  a  wave  model.  This  dissipation  is  incorporated  in  the  UCWP  model 
to  consume  wave  energy  corresponding  to  an  exponential  decay  rate.  The  non- 
linear wave-wave  interactions  are  included  in  the  UCWP  model  based  on  a  mod- 
ified JONS  WAP  spectrum  proposed  by  Donelaii  et  al.(l985).  The  migration  of 
peak  energy  frequency  caused  by  this  wave-wave  interaction  mechanism  is  deter- 
mined in  the  model  based  on  the  idealized  prognostic  equations  derived  by  Hassel- 
mann  et  al.(l976).  The  threshold  spectrum  which  indicates  a  minimal  condition  for 
the  occurrence  of  nonlinear  wave-wave  interactions  is  empirically  determined  based 
upon  the  field  data. 

The  UCWP  shallow  water  submodel  is  developed  with  the  inclusion  of  most 
of  the  known  significant  shallow  water  effects,  such  as  refraction,  shoaling,  perco- 
lation, bottom  friction,  and  wave  breaJcing  due  to  depth  limitation,  etc.  All  these 
shallow  water  effects  have  been  documented  in  the  past  (Chen  and  Wang,  1983; 
Shemdin  et  al.,  1978).  The  turbulent  dissipation  of  wave  energy  is  also  included  in 
the  model.  The  nonlinear  wave-wave  interaction  is,  however,  not  included  in  the 
model.  The  growth  of  waves  due  to  wind  input  is  evaluated  in  the  model  based  on 
a  TMA  spectrum,  as  suggested  by  Bouws  et  al.(1983).  The  five  parameter  TMA 
spectrum  has  been  correlated  to  wind  and  wave  conditions  by  Vincent  (1984). 

During  the  autumn  and  winter  seasons  extratropical  storms  and  hurricanes  fre- 
quently occur  in  the  Northwest  Atlantic  Ocean.  Waves  generated  by  the  storms  and 


6 

hurricanes  in  general  axe  the  major  cause  of  severe  coastal  erosion  along  the  south- 
east coast  of  the  United  States.  The  model's  performance  is,  thus,  tested  against 
such  severe  weather  conditions.  The  results  obtained  from  the  UCWP  deepwater 
submodel  axe  compared  with  measiired  wave  data  from  NOAA  buoys  deployed  in 
the  Northwest  Atlantic  Ocean.  The  results  obtained  from  the  UCWP  shallow  water 
submodel  axe  compared  with  measuied  wave  data  from  CDN  wave  gages. 

In  sumxnary,  the  detailed  research  included  in  this  study  can  be  briefly  described 
as  follows: 

(1)  A  niimber  of  standard  spectral  functions  and  angular  spreading  distribution 
functions  axe  reviewed.  These  axe  important  in  defining  an  appropriate  general 
directional  wave  spectrum. 

(2)  Two  existing  deepwater  wave  energy  spectral  models,  the  discrete  model 
by  the  Wave  Study  Information  (WIS)  group  of  Waterways  Experiment  Station 
(Resio  and  Tracy,  1983)  and  the  parametric  model  by  Hasselmann  et  al.(l973),  are 
rederived  to  highlight  the  basic  assumptions  and  fundamental  principles. 

(3)  Deepwater  wave  generation  mechainisms  are  reevaluated  based  upon  the 
currently  acquired  wave  data. 

(4)  Equations  that  are  either  theoretically  derived  or  empirically  proposed  as 
shallow  water  source  functions  with  respect  to  wave  generation,  propagation,  and 
dissipation,  are  examined  and  summarized.  These  are  later  utilized  to  build  a 
rigorous  shallow  water  model. 

(5)  A  new  discrete  model,  naxned  here  as  the  Unified  Coastal  Wave  Prediction 
(UCWP)  model,  is  developed  for  generating  wave  spectra  in  the  deep  ajnd  shallow 
water.  The  model  links  a  wind  field  model  to  a  combined  deepwater  and  shallow 
water  wave  prediction  model. 

(6)  Computer  algorithm  for  the  UCWP  model  is  developed  to  facilitate  numer- 
ical computations. 
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(7)  Results  from  the  deepwater  UCWP  model  are  compared  with  the  SMB 
model,  the  WIS  model,  and  Hasselmann's  (1976)  parametric  model  for  fetch-limited 
ajid  duration-limited  Ccises. 

(8)  Model  outputs  are  compajed  with  those  measured  by  deepwater  buoys  and 
shallow  water  wave  gages. 

(9)  Relative  importajice  of  individual  elements  affecting  the  deep  and  shallow 
water  spectral  transformations  is  analyzed. 

(10)  Different  values  of  bottom  friction  coefficient  are  tested  in  the  UCWP 
shallow  water  model  to  investigate  the  sensitivity  of  model  responses. 


CHAPTER  2 
REVIEW  OF  EXISTING  MODELS  AND  THEIR  PERFORMANCE 


2.1     Introduction 

The  approach  taken  in  the  present  study  for  ocean  wave  hindccisting  along  the 
Florida  coast  is  to  develop  a  numerical  model  based  on  state-of-the-art  deepwater 
aind  shallow  water  wave  prediction  techniques.  Currently,  deepwater  wave  models 
are  either  of  discrete  type  which  calculate  wave  energy  spectra  based  on  discretized 
wave  energy  equation  or  parametric  type  which  calculate  a  set  of  parameters  that 
in  turn  specifies  the  energy  spectra.  The  existing  shallow  water  wave  models  are 
invariably  of  discrete  type  since  the  mechanisms  involved  in  the  wave  transport  in 
the  shallow  water,  such  as  wave  refraction,  shoaling,  percolation,  wave  energy  dis- 
sipation, and  shallow  water  wave-wave  interactions,  are  too  complicated  to  permit 
parameterization. 

The  material  presented  in  this  chapter  includes  the  definitions  of  directional 
spectra  as  well  as  the  unidirectional  frequency  spectra  which  play  essential  roles 
in  wave  hindcast  models,  and  a  presentation  of  the  typical  theories  and  methods 
applied  to  the  ciirrent  deepwater  and  shallow  water  wave  models.  Two  widely  known 
deepwater  wave  models,  the  parametric  spectral  model  by  Hasselmann  et  al.(1973) 
and  the  discrete  spectral  model  by  the  Wave  Information  Study  (WIS)  group  of 
the  Waterways  Experiment  Station  (WES),  are  reviewed  as  a  demonstration  of 
spectral  methods.  It  is  shown  that  the  deepwater  parametric  spectral  model  in 
general  is  meant  for  wind  waves  only  and  cannot  account  for  the  swells  while  the 
deepwater  discrete  spectral  model  by  WIS  has  no  difficulty  in  tracing  both.  It  is 
further  shown  that  the  WIS  spectral  model  has  the  same  ability  as  the  parametric 
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model  to  include  the  nonlinear  wave-wave  interactions  in  computing  wave  energy 
transformation  in  the  deep  water;  this  is  accomplished  using  the  formula  which 
estimates  the  nonlineair  wave  energy  tremsformation  rate  in  the  frequency  domain 
derived  from  the  prognostic  equations  of  pajametric  method. 

The  shallow  water  wave  model  is  reviewed  following  the  works  done  by  Chen 
and  Wang  (1983),  Collins  (1972),  and  Shemdin  et  al.(1978).  Among  the  major 
contributors  of  shallow  water  effects,  wave  refraction  and  shoaling  axe  included  in  the 
advection  terms,  and  others  axe  acting  as  a  part  of  the  forcing  functions  governing 
the  generation,  dissipation,  and  nonlinecir  trajisfer  of  wave  energy  in  shallow  water 
region. 

The  model  performance  is  only  shown  here  for  the  deepwater  model  of  WIS. 
The  shallow  water  model  which  requires  the  output  of  deepwater  model  to  specify 
the  boimdary  condition  is  not  tested  until  the  latter  is  modified.  Its  performance 
will  be  discussed  in  Chapter  4. 

2.2     Spectral  Definitions 

A  directional  wave  spectrum  is  essentially  a  representation  of  the  wave  num- 
ber spectrum  of  water  surface  waves  that  follow  the  linear  wave  theory  dispersion 
relationship  between  wave  number  and  wave  frequency.  By  linear  superposition  of 
wave  components,  an  irregular  water  surface  can  be  approximated  by  taking  the 
stochastic  integral  over  the  wave  number  plane,  s,  which  is  taken  here  as  identical 
to  the  still  water  surface  plane  (Thomson,  1972). 

T]{si,t)==Re[ll  da(^)e'(* ''-"')]  =  ff  -[da(it)e'(*"*'-^')  +  da*()t)e-'(* '"'-"*']      (2.1) 

where  r)  is  water  surface  elevation,  t  is  time,  i*i  =  {xi,yi)  is  the  horizontal  location 
of  interest  in  rectangular  coordinates,  a  is  radian  frequency,  k  =  {k-cos9,  k-sinO)  is 
a  vector  wave  number  with  magnitude  A;  =|^|  and  direction  9,  da{k)  indicates  an 
infinitesimal  wave  amplitude  (complex  number)  which  is  the  contribution  of  waves 
in  a  wave  number  interval  dk,  an  asterisk  appearing  in  superscript  indicates  the 
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complex  conjugate,  the  operator  Re[  ]  indicates  the  real  part  of  a  complex  number. 
The  wave  amplitude  component  da{k)  can  be  transformed  into  da{cr,9)  since  k  can 
be  expressed  by  a  direction  9  and  a  frequency  a.  To  a  first  approximation,  the  wave 
frequency  satisfies  the  well-known  linear  wave  theory  dispersion  equation 

a^  =  gk-tsLiihikh),     k  =\k\  (2.2) 

where  g  is  the  gravitational  constant  and  h  is  the  water  depth.  A  small  contribution 
of  the  water  surface  waves  representing  a  specific  vector  wave  number  is 

*  dr7(si,0  =  -[da(fc)e'(*^'"'-''*)  +  da*(fc)e-*(*-'"'-''*)]  (2.3) 

Thus,  the  mean  square  value  of  dry  can  be  evaluated,  which  has  a  form 

[dTi{sut)\^    =    -  llm  ^r^\{da{k)]h'^^'^''-°'^+2da{k)da'{k) 

4  T—oo  1  J-T/2 

+  [da*(^)]'e-'2(*'"'-<")}(it  =  -da(^)da'(fc) 
=    l\da{k)\^=-\da{a,9)\^^E{a,e)dade^dE{(j,e)       (2.4) 

where  T  is  the  time  duration  of  interest,  ein  overhead  baj  indicates  the  operation  of 
arithmatic  mean;  E{a,  9)  is  the  directional  frequency  spectrum  which  is  defined  as 

E{a,  9)  =  [drj{si,t)]yd<7d9  =  -  | da{a,  9)  \ydad9  (2.5) 

According  to  linear  wave  theory,  the  wave  energy  per  unit  water  surface  area  per 
imit  water  weight  is  equal  to  half  the  squared  value  of  wave  amplitude  (Ippen,  1966). 
Therefore,  Eq.(2.4)  specifies  a  quantity  proportional  to  wave  energy  associated  with 
a  particular  frequency  range,  dcr,  cind  directional  interval,  d^.  However,  Eq.(2.4) 
may  be  regarded  as  a  contribution  from  an  infinitesimal  amplitude  wave  with  ampli- 
tude da{a^9)  at  the  mid  da  and  d9.  Thus,  E{a,9)  is  interpreted  as  the  distribution 
of  wave  energy  per  unit  water  surface  area  per  unit  water  weight  per  unit  frequency 
per  unit  direction.  The  corresponding  one-dimensional  frequency  spectrum  £'(cr)  is 
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defined  by  integrating  E{a,0)  over  all  6. 

E{a)  =  r  E{a,6)d9  (2.6) 

J—ir 

The  maxginal  directional  spectrum  is  defined  as 

/•oo 

E[d)=j     E{a,9)da  (2.7) 

The  total  wave  energy  per  unit  water  surface  area  per  unit  water  weight  is  defined 
by  integrating  E{a)  over  all  a,  or  integrating  E{d)  over  all  9,  or  integrating  E{a,9) 
over  all  a  and  9.  That  is, 

/•cx)  fir  rtr      roo 

Etotai=        E{a)da=        E{9)d9  =         /     E{a,9)d(Tdd  (2.8) 

•'0  J -r  J  -X  Jo 

Since  the  spectral  density  defined  in  Eq.(2.5)  or  Eq.(2.6)  has  dimensions  of  rate 
of  change  of  wave  energy  density,  the  name  'power  spectral  density'  is  often  used. 
However,  an  alternate  name  'energy  spectral  density'  is  more  commonly  used  in 
coastal  engineering. 

2.3     Spectral  Shape  Function  Models 

In  the  study  of  wave  hindcast  through  spectral  models,  the  standard  spectral 
shape  functions  are  served  as  fundamental  structure  unit  in  building  the  wave  en- 
ergy spectral  states  in  the  sea.  Three  frequency  spectral  shape  functions  are  well- 
accepted  in  the  past:  (i)  Pierson-Moskowitz  (P-M)  spectrum,  (ii)  Kitaigorodskii 
spectrum,  and  (iii)  JONS  WAP  spectrum.  The  Pierson-Moskowitz  spectrum  (1964), 
EpMif),  is  given  by 

^''^(^^  =  (^-P(-0-74(/o//)^]  =  J^s^M-lifm/m  (2.9) 

where  /  =  cr/27r  is  the  wave  frequency  with  units  l/time,  Op  is  known  as  Phillips' 
constant,  /o  =  g/27r(7i9.5  is  the  frequency  of  the  deepwater  waves  whose  phase 
speed  is  identical  to  the  wind  speed  Uig,s  at  a  height  19.5  m  above  sea  surface, 
and  /„  —  0.88/0  is  the  wave  frequency  of  the  spectral  peak.  More  often,  the  name 
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'spectral  peak  frequency'  is  referred  to  the  frequency  of  spectral  peaJc,  /^,  and  the 

wave  period  corresponding  to  this  frequency  is  known  as  'modal  period'. 

The  Kitaigorodskii  spectrum(1962),  £'^(/),  is  given  by 

r  ag'f-'{2n)-*  exp[l  -  {fj /)'],    iU  <  fm 
EkU)  =  (2.10) 

i  ag'f-'{27r)-\  li  f  >  U 

which  contains  two  free  parameters,  the  energy  scale  factor  a  and  the  spectral 
peaJc  frequency  fm-  The  JONSWAP  spectrum  is  derived  from  the  Joint  North  Sea 
Wave  Project  (Hasselmcinn  et  al.  1973)  and  constitutes  a  modification  to  the  P-M 
spectnira  to  provide  for  a  much  more  sharply  peaiked  spectrum.  The  JONSWAP 
spectriim,  Ej{f),  is  given  by 

aa^  5  exp[--(^^-^^)2] 

with 

_    /    <7a,       li  f  <   fm 
•^  ~    1    ^6,       if  /   >  /. 

The  function  contains  five  free  parameters,  two  scale  parameters  a  and  /„,  and 
three  shape  pairameters  7,  aa  and  Cb.  The  parameter  7  is  the  peak-enhancement 
factor  which  is  the  ratio  of  the  peak  value  of  the  spectrum  to  the  peak  value  of  the 
corresponding  P-M  spectrum  with  the  same  values  of  a  and  fm-  The  parameters 
Ca  and  a),  yield  the  scales  of  left  and  right  peak  widths,  respectively. 

Figure  2.1  shows  a  comparison  of  spectral  shapes,  normalized  by  the  spec- 
tral peak  energy  of  the  P-M  spectrum,  of  the  P-M  spectrum,  the  Kitaigorodskii 
spectrum,  and  the  mean  JONSWAP  spectrum,  which  has  its  spectral  parameters 
7,  Oa,  and  (Ti  equal  to  3.3,  0.07,  and  0.09,  respectively.  It  is  seen  that  the  mean 
JONSWAP  spectrum  and  the  Kitaigorodskii  spectrum  have  similar  shape  and  they 
are  much  more  sharply  peaked  than  the  P-M  spectrum. 

All  three  spectra  shown  above  indicate  an  'equilibrium  range'  on  the  rear  face 
of  the  spectral  shape  with  the  associated  spectral  energy  densities  being  inversely 
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Figure  2.1:  Comparison  of  nondimensionalized  shapes  of  the  P-M  spectrum,  the 
Kitaigorodskii  spectrum,  the  JONSWAP  spectrum  ajid  the  modified  JONSWAP 
spectrum. 
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proportional  to  frequency  to  the  fifth  power.  However,  more  recent  studies  made  by 

Donelan  et  al.(l985)  and  Battjes  et  al.(1987)  using  both  field  and  laboratory  data 

showed  that  an  /"■*  frequency  dependence  at  the  tail  part  of  the  wave  spectrum 

may  be  more  adequate  for  the  wind-generated  waves  than  the  f~^  spectral  shape. 

Based  upon  the  assvmaption  of  the  /"■*  dependence  on  the  rear  face  of  frequency 

spectrum,  Donelaji  et  al.(1985)  suggested  a  modified  JONS  WAP  spectrum 

!        /  J        -/"i  ^2l 
nn^  exp ( ;: — ) 

EjmU)  =  j^^^^exv[-{fm/mi  2    a.,/^  (2.12) 

with  all  its  free  parameters  a,  /„i,  7,  and  crab  defined  the  same  as  in  JONSWAP 
spectrum.  This  modified  JONSWAP  spectrum  is  immediately  applicable  to  any 
discrete  spectral  method  provided  fm  is  known  or  defined  in  advance.  This  is 
because  all  the  empirical  relationships  between  each  of  the  free  parameters  a,  -7, 
and  Oab  were  established,  based  on  field  data,  by  Donelan  et  al.(1985)  as  functions  of 
fo/ fm-  A  modified  JONSWAP  spectrum  with  its  parameters  7  and  aaj  being  equal 
to  the  mean  JONSWAP  ones  is  also  shown  in  Fig.  2.1.  This  spectrum  hcis  wider 
width  and  higher  peak  than  the  mean  JONSWAP  spectrum,  and  it  contains  more 
energy  in  the  high  frequency  components  than  the  mean  JONSWAP  spectrum,  the 
P-M  spectrum  and  the  Kitaigorodskii  spectrum. 

The  directional  dependence  of  the  wave  spectrum  is  usually  modelled  by  a  sym- 
metrical directional  distribution  fiinction  which  is  also  known  as  the  'angular  distri- 
bution function'  or  'spreading  function'.  It  is  often  convenient  to  define  a  normalized 
spreading  function  as  the  ratio  of  (Cartwright,  1963) 

H.(/.»)  =  ^  (2.13) 

where  E{f,6)  =  Hs{f,0)E{f)  indicates  the  directional  spectrum,  which  has  its 
wave  energy  densities  distributed  in  the  {f,0)  domains  according  to  Hs{f,0). 

A  number  of  functional  forms  of  Hs{f,0)  have  been  suggested  but  all  with 
similar  shape,  i.e.,  they  are  bell-shaped  and  symmetric  about  a  central  direction.  It 
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Figure  2.2:  Plots  of  Hs{f,9)  versus  0  -  ^o- 

is  usually  difficult  to  choose  one  in  preference  to  another  in  practical  applications. 
A  commonly  used  function,  which  is  proportional  to  a  power  cosine  function,  is 
given  here  (Sand,  1979). 

TV   r(25  +  l)  ^2.14) 

0,  elsewhere 


Hs{f,e)  = 


where  5  =  S{f)  is  called  the  spreading  parameter,  which  is  always  positive  and  a 
function  of  frequency,  ^o  =  ^o(/)  is  the  central  angle  of  the  symmetric  spreading 
directions.  Figure  2.2  shows  several  distributions  of  Hs{f,0)  with  the  spreading 
parameter  ranging  from  1  to  20.  It  is  shown  that  Hs{f,0)  associated  with  higher 
spreading  parameter  can  have  very  sharply  peaked,  narrow-band  directional  distri- 
bution. When  this  Hs{f,9)  is  applied  in  a  wave  prediction  model,  the  spreading 
parameter  S  is  commonly  defined  equal  to  1  and  2,  respectively,  for  all  frequency 
components  of  locally  generated  and  nonlinearly  transfered  energies.     It  can  be 
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seen  in  Fig.  2.2  that  both  spreading  functions,  Hs{f,0),  corresponding  to  5  =  1 
and  5  =  2,  respectively,  have  relatively  wide  distributions  over  a  180°  directional 
domain. 

2.4     Deepwater  Pajaimetric  Model 

In  the  absence  of  current  in  the  ocean,  the  general  equation  governing  the  trans- 
port of  wave  energy  is  given  by  (Hasselmajin  et  al.,  1973) 

,_EU^  ^  em^  ^  C.if,s-f-^^l^  =  G(/.M- 0  (2.15) 

where  E{f,0,s,t)  is  the  energy  density  of  the  spectral  element  with  reference  to 
still  water  surface  plame  at  horizontal  location  s  and  time  t,  with  frequency  /  and 
direction  of  propagation  9,  i  is  the  path  of  this  spectral  element  in  which  d£  = 
Id^le'*,  Cg{f,s)  =  da/dk\s-  =  d[2TTf)ldk\s  is  the  group  velocity  associated  with  this 
spectral  element,  k  is  the  magnitude  of  the  associated  wave  nimiber,  and  G{f,6,s,t) 
is  the  net  source  function  representing  the  rate  of  energy  transfer  into  or  out  of  this 
element.  This  equation  is  commonly  referred  to  as  'wave  energy  balance  equation', 
or  'wave  energy  transport  equation',  or  simply  'transport  equation'.  In  the  text  to 
follow,  the  notations  of  E{f,  6)  and  G{f,  9)  are  often  used  for  simplicity. 

The  magnitude  of  wave  group  velocity  Cg{f,s)  can  be  derived  from  the  disper- 
sion relation  given  in  Eq.(2.2)  and  expressed  as 

C.(/,.-)=n.C  =  (i.-^)C  (2.16) 

where  C  =  a/k  =  2nf/k  is  phase  speed  of  waves  with  frequency  equal  to  /,  h  is 
water  depth  at  the  location  5,  and  n  is  the  ratio  of  wave  group  velocity  to  wave 
phase  velocity.  It  is  noted  that  both  phase  velocity  C  aind  wave  number  k  aire 
implicit  functions  of  /.  The  presence  of  the  source  function  G{f,9,s,t)  is  thought 
to  include  all  the  possible  physical  terms  which  are  responsible  for  the  local  growth 
auid  decay  of  water  surface  waves. 
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The  deepwater  wave  is  usually  defined  in  coastal  engineering  by 

tanh(A;/i)  >  0.9962,     or    {kh)  >  it  (2.17) 

Under  the  above  condition,  the  wave  properties  are  hajdly  influenced  by  the  bottom. 
Therefore,  the  shallow  water  effects,  such  as  shoaling,  refraction,  wave  dissipation 
due  to  bottom  friction,  etc.,  are  no  longer  important  in  the  deepwater  wave  hindcast. 
It  is  noted  that  in  deep  water,  the  dispersion  relation  is  reduced  to 

a'  =  {2nfY  =  gk 

emd  the  deepwater  wave  group  velocity,  Cgo{f),  can  be  expressed  as 


2  2a       47r/ 

which  is  not  dependent  on  water  depth  and,  hence,  is  not  dependent  on  location. 

The  parametric  spectral  model  was  originally  developed  by  Hasselmann  et  al. 
at  1973.  The  model  was  later  applied  by  many  investigators  without  further  mod- 
ifications (Weare  and  Worthington,  1978;  Gunther  et  al.,  1979).  The  presentation 
of  the  parametric  model  in  this  section  is  mainly  based  on  the  work  of  Hassel- 
mann et  al.(1973). 

To  derive  the  equations  solving  the  spectral  parameters,  the  directional  wave 
spectrum  used  in  the  transport  equation,  Eq.(2.14),  is  assumed  to  have  the  following 
expression: 


Ej{f,d)  =  Ej{f)H,{f,9) 


2 

Ej{f)-cos'^{9  -60),        if  1^ -^o|<7r/2 

^  (2.18) 


0,  elsewhere 

where  Hi{f,8),  the  spreading  function,  is  a  special  case  of  Eq.(2.14)  with  S  equal 
to  1,  and  ^0  coincides  with  the  downwind  direction.  By  first  substituting  Eq.(2.18) 
into  the  transport  equation,  Eq.(2.15),  and  then  intergrating  the  transport  equation 
over  the  direction  domain,  0  <0<  2n,  the  direction-independent  transport  equation 
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is  obtained  as 


^  =  ^.^(/,^  =  .(/.  (...) 


with  the  directional  averaged  group  velocity  Cgo{f)  equal  to 

LV„{f)EAMM  _    8     5 


CM  =  -'  '-{;{'   ■       =  ~  (2.20) 

ajid  the  directional  integrated  source  function  G[f)  equal  to 

G{f)  =  l'G{f,e)dd  (2.21) 

Note  that  symbols  s  and  t,  which  stajid  for  location  and  time,  respectively,  are 
not  shown  in  Eqs.(2.19),  (2.20),  zmd  (2.21),  and  in  the  following  derivations  for 
simplicity. 

The  source  function  G{f)  is  usually  represented  by  three  quantites,  the  energy 
input  from  the  atmosphere,  Gin{f),  the  dissipation  of  energy,  Gdi,{f),  and  the 
nonlinear  energy  transfer,  G^(/).  That  is, 

G{f)  =  G.„(/)  +  G,,,(/)  +  Gr^tif)  (2.22) 

Since  the  spectral  shape  is  assumed  to  be  consistent  with  the  JONSWAP  spectrum 
in  the  parametric  model  and  is  presumed  in  a  state  of  equilibrium  range,  the  dis- 
sipation Gdi,{f)  is  in  general  dropped  from  Eq.(2.22).  The  rest  two  source  terms 
Gin{f)  asid  Gni{f)  are  thought  to  be  equally  important  for  the  growth  and  transfer 
of  wave  energy  in  the  ocean.  It  was  demonstrated  by  Resio  and  Vincent  (1979)  in 
numerical  model  studies  that  the  duration  required  to  achieve  fetch-limited  condi- 
tion in  the  case  dominated  by  atmospheric  input  is  of  the  same  order  eis  the  case 
dominated  by  the  wave-wave  interactions. 

The  growth  rate  G,n(/)  is  conmionly  approximated  in  terms  of  a  linezir  func- 
tion of  one-dimensional  frequency  spectrum  (Snyder  and  Cox,  1966).  However,  for 
convenience  the  quantity  G,„(/)  is  often  assumed  to  be  proportional  to  E{f)  in  the 
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parametric  model.  For  instance,  the  energy  input  source  term,  utilized  by  Hassel- 
mann  et  al.(1973)  in  their  wave  model  is  linearly  proportional  to  one-dimensional 
frequency  spectrum  and  has  the  following  expression: 

G,M)  =  5sf^Ej{f)  (2.23) 

where  sss  0.0012  is  the  ratio  of  air  density  to  sea  water  density,  and  Uw  is  the  wind 
speed  at  10  m  elevation  above  sea  surface. 

The  general  form  of  the  nonlinear  transfer  source  Gni{f)  for  a  single  peaked 
spectrum  is  known  to  have  a  positive  lobe  appearing  on  the  frontface  of  the  spec- 
trum, a  negative  second  lobe  appearing  on  the  rearface  of  the  spectrum  closer  to  the 
spectral  peak,  and  a  positive  third  lobe  appearing  on  the  reaxface  of  the  spectrum 
away  from  the  spectral  peak  (see  Fig.  2.3).  In  wave  hindcast,  only  the  first  two 
lobes  are  important  since  the  energy  transfer  into  the  third  lobe  is  at  such  frequen- 
cies that  viscous  and  turbulent  dissipations  are  cissumed  to  keep  the  local  energy 
components  small. 

The  nonlinear  transfer  function  deduced  by  Hasselmann  et  al.(1973)  has  the 
expression 

2     3 

Gniif)  =  Di  ^-^rPif/fm)  (2.24) 

Jm 

where  Z?i  is  a  proportional  constajit  whose  magnitude  is  on  the  order  of  10"^, 

and  ^{f/fm)  is  a  non-dimensional  function  depending  on  the  spectral  shape.    By 

neglecting  the  energy  transfer  in  the  third  lobe,  where  waves  of  high  frequencies  are 

more  likely  governed  by  the  equilibrium  f~^  law,  Eq.(2.24)  can  be  approximated  by 

the  expression 

2.2a'fmEj{f),        \i  f  <fm 
Gniif)  =  (2-25) 

-2.2a'UEj{f),       if  f>U 

With  the  source  function  given,  the  direction-independent  transport  equation, 
Eq.(2.19),  can  be  then  projected  onto  the  JONS  WAP  parameter  space  to  yield 
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Figure  2.3:  Nonlinear  transfer  of  energy  (Hasselmann  et  al.,  1973). 
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(Hasselmajin  et  ah,  1973) 

dbi       ^    dbi       ^  .    .  ,        , 

-^  +  Ai^  =  Gi,        t,j  =  1,2,3,4,5  (2.26) 

where  6i  =  /m>  &2  =  a,  63  =  7,  b^  =  (Ja,  and  65  =  aj.  The  equations  given  in 
Eq.(2.26)  eire  known  as  the  'prognostic  equations'  in  the  parametric  model.  The 
Dij  matrix  is  a  generalized  propagation  velocity,  ajid  G,-  are  the  generalized  source 
functions.  A  list  of  the  Dij  matrix  and  G,-  functions  so  obtained,  with  Gdi,[f)  equal 
to  zero,  is  given  in  Appendix  A. 

Thus,  by  solving  Eqs.(2.26)  simultcineoiisly  with  given  wind  information,  the  five 
free  JONSWAP  parameters  can  be  obtained  and  the  resulting  JONSWAP  spectrum 
can  be  then  expanded  according  to  the  Hi{f,9)  spreading  function  to  generate  a 
directional  distribution  of  wave  energy  densities. 

The  nature  of  the  parametric  approach  is  that  only  single-peaked  spectra  are 
allowed  in  the  hindcast,  and  as  a  result  temporal  and  spatial  variations  of  swells  can 
not  be  included  in  the  model.  To  include  the  swell  portion  of  the  wave  energy,  it  is 
necessary  to  incorporate  a  discrete  spectral  method  to  accommodate  the  transport 
of  swells.  Unlike  the  parametric  method,  which  resorts  to  solving  the  free  spectral 
parameters  of  a  predesignated  spectral  function,  the  discrete  spectral  method  tracks 
the  spectral  densities  of  both  the  preexisting  and  currently  generated  waves  in  the 
wave  field  of  interest.  Accordingly,  the  spectra  so  obtained  could  have  arbitrary 
shapes  including  multiple  peaks.  The  details  of  the  discrete  spectral  method  are 
given  in  the  next  section. 

2.5     Deepwater  Discrete  Model 

The  development  of  deepwater  discrete  spectral  models  is  again  bjised  on  the 
wave  energy  transport  equation  given  in  Eq.{2.15).  The  method  utilized  in  discrete 
spectral  models  generates  directly  the  directional  distributions  of  wave  energy  ac- 
cording to  the  local  atmosphere  energy  input  and  then  transfers  them  advectively  in 
the  horizontal  domain  according  to  the  governing  equation.  Since  the  fundamental 
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structure  of  a  discrete  spectral  model  does  not  change  from  one  model  to  another, 
the  solution  procedure  used  in  a  discrete  spectral  model  will  be  described  in  this 
section  through  the  presentation  of  the  discrete  model  developed  by  Wave  Informa- 
tion Study  (WIS)  group  at  the  Waterways  Experiment  Station  (WES),  U.S.  Army 
Corps  of  Engineers. 

The  WIS  model  was  originally  created  to  hindcast  the  deepwater  waves  in  the 
Atlantic  Ocean  between  1956  amd  1975.  The  documentation  of  the  model  can  be 
found  in  the  series  of  WIS  reports  of  Vol.  8  (Jensen,  1983a)  and  12  (Resio  and 
Tracy,  1983).  The  version  of  the  WIS  model  presented  here  is  the  one  modified  by 
the  same  WIS  group  in  1984  by  studying  the  deepwater  wave  hindcast  in  the  Pacific 
Ocean. 

The  WIS  model  solves  Eq.(2.15)  numerically  based  on  the  following  discretized 
form: 

E^^'ifJ)  =  E\f,e)  -  C,,(/)^:^M  .  Af  +  G{f,d)  ■  At  (2.27) 

where  the  superscript  indicates  the  time  level.  The  second  term  on  the  RHS  is 
energy  advection  and  the  third  term  is  the  net  energy  source,  and  these  two  terms 
axe  uncoupled. 

The  advective  transfer  of  wave  energy  is  computed  in  the  model  using  the 
method  developed  by  Hassehnann  et  al.(1968).  This  method  divides  the  entire 
direction  domain,  0  <  ^  <  27r,  uniformly  into  sixteen  directional  bands,  each  hav- 
ing an  angular  spaji  of  27r/16  radian,  and  wave  energy  contained  in  each  band  is 
advected  independently  without  coupling  with  others.  To  express  the  method  in  a 
numerical  scheme,  the  advective  energy  E{f,0)ad  may  be  written  explicitly  as 

E^^'ifJ)  =  E\f,9)  -  C,4f)^^^^  .  At 

The  source  term  is  further  divided  into  three  components: 

G{f,9)  .  At  =  [G^:HfJ)  +  G',„{f,0)  +  G^M^O)]  .  At  (2.28) 


23 


where  G,„,  Gdis,  and  Gru  represent  atmospheric  energy  input,  wave  energy  dissi- 
pation and  nonlineax  wave-wave  interaxition,  respectively.  The  G,„  term,  which  is 
computed  at  the  current  time  level,  is  further  split  into  two  parts:  (i)  G,„_i(/,  ^), 
the  rate  of  growth  of  the  prevailing  waves  and  (ii)  G,n,2(/?  ^)>  the  rate  of  generation 
of  new  wind  waves  in  the  downwind  direction. 

In  WIS  model,  the  two  atmospheric  input  terms  are  formulated  as  follows: 

g'a'       8 


G.n.iAi  = 


4.4  Atf„ 


0, 


(27r)V^  37r 


cos 


'{0-9),        if/>/„ 


elsewhere 


with  the  mean  propagation  direction,  9,  defined  as 

-_        _,Jfsm9E{f,9)dfd9 
'-"■^     ^IJcos9E{f,9)dfd9- 

where  the  scale  parameters  a  and  /^  are  to  be  determined  later. 


} 


(2.29) 


(2.30) 


Gin,  2  At 


-^^exp[l  -  {l^y]^cos\9  -  9^),     if/.  >  0.25  (Hz),  /  <  U, 
(27r)V^  /       TT 

f^   >  0.2  (Hz) 

if/,  >  0.25  (Hz),/ >/^, 


ag 


2      2 


(27r)V'7r 
0, 


cos^9-6^), 


U   >  0.2  (Hz) 
elsewhere 


(2.31) 


where  0^  denotes  downwind  direction  and  /,  denotes  a  peak  energy  saturation 
frequency  which  is  not  necessarily  equal  to  the  peak  spectral  frequency  and  is  defined 
as 


E{f,)  =  max  {E{f)  >  E.{f)} 
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Ea{f)  is  the  one-dimensional  saturation  spectrum  defined  as 


""^^  if  /  <  /o  -  g/2irUw 


Esif)  = 


(27r)V^   ' 

J^pi^  -  ifo/fmV]  +  Ep^{f){fol U)\    elsewhere 


(2.32) 

where  /o  is  the  wind-wave  resonajit  frequency  to  which  the  corresponding  wave 
phase  speed  is  identical  to  the  wind  speed  at  10  m  elevation,  C/io,  and  Epf^  denotes 
the  P-M  spectrum. 

There  are  two  scale  parameters  a  and  fm  to  be  determined  first  at  time  level 
k+1.  The  scale  parajneter  f^.  is  obtained  by  solving  the  prognostic  equations  de- 
rived in  the  (a,/^)  plajie  a^  described  by  Hasselmann  et  al.(1976)  in  their  para- 
metric model.  By  allowing  variations  of  only  two  free  parameters  a  and  /^  in  the 
JONSWAP  parametric  spax;e,  the  a—  ajid  fm—  component  prognostic  equations  can 
be  derived  as 

^  +  0.47  C;;(/)^  +  O-^CiM)^^  -  0.005  i^-^r^'afm  -  5  a' f^    (2.33) 


and 


^^"^  +C',.(/)^  -0.07C,„(/)^^  =  -0.54aV^  (2.34) 


If  the  changing  wind  speed  can  be  expressed  as  a  power  law  of  time, 


U,o 

-</ 

with  Uq 

a  reference 

wind  speed. 

fm  < 

ind  a  can  be  sc 

>lved  from 

Eqs.(2.33) 

and 

(2.34) 

as 

g 

'Jo 

-1) 

(2.35) 

and 

a  - 

-  5f  ■'"'^^°)2/3 

(2.36) 
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with,  constants  A  and  B  given  by  (Hasselmann  et  al.,  1976) 

A  =  16.8  (1  +  1.51  q)^^\      B  =  0.031  (-^-^ -Y^^ 

1  +  1.51  q 

To  obtain  fm  at  time  level  k+1,  the  value  at  t  in  Eq.(2.35)  is  replaced  by  t  +  At, 
which  results  in 

^SlllSl  =  A(fl  +  ^)(3/7)(,-i)  (2.37) 

Eliminating  gtlUo  from  Eqs.(2.35)  and  (2.37)  yields 

^/^■"^  =  {[fl)^'"^^'-'^  +  ^(^)-(V3)(,-l)}(3/7)(,-l)  (2.38) 

U-iQ  Uq     Ag 

For  constajit  wind  speed  within  the  duration  At,  the  above  equation  reduces  to 

f'J'  =  iifi)-'''  +  a^ii^y^'At}-'/'  (2.39) 

where  a^  is  the  nonlineax  wave  energy  transport  coefficient  defined  by  Hcissel- 
mann  et  al.(1976)  ajtid  suggested  as  equal  to  0.00138.  In  WIS  model,  the  f^^ 
value  is  replaced  by  /o  if  the  spectral  peaJc  frequency  is  less  thaji  /q.  Also,  instead 
of  using  Eq.(2.36)  to  compute  a,  WIS  model  employs  the  following  relationship 
suggested  by  Resio  and  Vincent  (1977): 

a  =  0.0489  (^-^)-°-2  (2.40) 

where  U,  is  the  surface  frictional  velocity  obtained  by  assuming  logarithmic  over 
water  wind  profile 

^  =  -ln-  (2.41) 

U,       0.4       20  ^        ' 

where  U^  is  wind  speed  at  a  height  z  above  the  sea  surface,  and  zq  is  the  equivalent 
roughness  at  water  surface  given  by 

0.1525  Uj  ,  X  ,        V 

Zq  =  —- —  +  0.0144—^  -  0.00371,      (in  cgs  units)  (2.42) 

u*  g 
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Although  Eq.(2.36)  is  not  used  in  the  WIS  model,  it  can  be  linked  to  an 
empirical  relationship  {g^Etotai/U*o){f„iUio/g)*/a=consta.nt  suggested  by  Hassel- 
mainn  et  al.(1976)  to  yield  a  power  law  a  ~  {g^Etotcu/UiQ)'^'^  which  is  similar  to 
Eq.(2.40). 

The  transfer  of  wave  energy  due  to  a  nonlinear  wave-wave  interaction  source 
Gruifi^)  is  formulated  based  upon  Kitaigorodskii  spectrum  (1962)  and  a  spreading 
fiinction  with  the  central  angle  parallel  to  the  mean  wave  direction  d  (Resio,  1981). 


Gni^t  =  { 


0.0023  At^rPif/fm)  ^  cos^(5  -  0),        if  |^  -  ^|  <  7r/2 

J  171 


(2.43) 


0,  elsewhere 

where 

f    exp[l  -  ifm/m,  If    f<fm 

[    -[fmlf)\  if    /  >   /m 

The  dissipation  includes  two  parts:  the  energy  loss  of  wind  waves  and  the  energy 
loss  of  swell.  The  energy  loss  of  wind  waves  is  accounted  for  by  restricting  the  growth 
to  saturated  spectral  condition  with  f~^  law  on  the  downslope  side.  The  loss  of  wave 
energy  in  swell  is  adopted  from  Resio  (1981),  in  terms  of  one-dimensional  frequency 
spectrum. 

G,.-,Ai  =  ^^""^  3  (2.44) 

0,  elsewhere 

where  f^e  indicates  peak  spectral  frequency  of  swell. 

2.6     Shallow  Water  Model 

As  deepwater  waves  propagate  into  area  of  shallow  water,  they  will  sense  the 
bottom  and  wave  energy  will  be  redistributed  asid  dissipated  in  the  direction  and 
frequency  domains  due  to  refraction,  shoaling,  bottom  friction,  percolation,  etc. 
Among  these  shallow  water  effects,  only  refraction  and  shoaling  can  be  handled 
fairly  readily  in  many  situations,  others  require  further  theoretical  and  experimental 
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studies.   Shallow  water  waves  can  also  grow  if  winds  are  present  and  large  enough 

to  add  energy  into  waves  as  a  result  of  air-sea  interaction. 

In  shallow  water  the  wave  energy  transport  equation  can  be  expressed  as  follows: 
(Chen  and  Wang,  1983) 

dE{f,9)    ,   dE{f,9)C,if,e)cose  ,   dEjf, 9)C,{f, 9)  sme 

—ar-  + di "^ d^ =  ^^^'  ^^      ^^-''^ 

where  9  is  defined  counterclockwise  from  the  positive  x  axis.  Here,  unlike  the 
deepwater  case,  both  Cg  ajid  9  are  no  longer  independent  at  position  (actually 
water  depth).  In  order  to  track  energy  trajisport,  changes  in  wave  direction  due  to 
wave  refraction  must  be  considered  simultaneously.  This  requires  solving  additional 
equation  governing  the  chajige  of  wave  number  vector.  According  to  the  refraction 
law,  the  curl  of  wave  number  is  identical  to  zero,  or 

„      -.       (dkcos9)       (dksm9) 

VxA;=^— --^ ^  =  0  2.46 

oy  ox 

Therefore,  if  the  sotirce  function  G{f,  9)  is  known  in  advance,  the  directional 
spectrum  E{f,9)  in  shallow  water  can  be  obtained  by  solving  both  Eqs.(2.45)  and 
(2.46).  It  should  be  noted  that  neither  current  nor  tide  is  considered  to  influence 
the  wave  spectral  transformation  governed  by  Eqs.(2.45)  and  (2.46).  Although  it  is 
possible  to  include  the  current  and  tide  effects  in  solving  the  wave  spectra  in  the 
shallow  water  (Chen  ajid  Wang,  1983),  the  influences  of  currents  and  tides  to  the 
wave  spectral  transformation  in  the  shallow  water  will  be  exclusive  in  the  present 
study  for  the  following  two  reasons.  First,  the  numerical  model  which  includes  the 
current  and  tide  effects  in  estimating  shallow  water  wave  spectra  is  in  general  dealing 
with  a  large  matrix  of  equations  such  that  execution  of  the  model  is  not  efficient 
in  the  aspect  of  computational  time.  Second,  the  insufficient  information  about 
interactions  among  waves,  currents,  and  tides,  according  to  the  present  knowledge, 
may  cause  problem  in  defining  the  source  terms  associated  with  both  current  and 
tide  influences. 
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The  source  function  G{f,9)  may  still  be  represented  by  the  summation  of  three 
contributions,  namely,  generation,  dissipation,  and  nonlinear  wave-wave  interaction, 
although  the  formulais  are  expected  to  be  different  from  the  deepwater  ones  due  to 
the  influence  of  the  bottom. 

G{f,0)  =  G,4f,0)  +  G,u{f,e)  +  Gr,t{f,9) 

Explicit  formulas  for  individual  source  terms  are  well  documented  in  the  pcist  (Chen 
and  Wajig,  1983;  Collins,  1972;  Shemdin  et  al.,  1978)  and  are  briefly  described  in 
the  remainings  of  this  section. 

The  generation  term  is  usually  defined  by  a  combined  Phillips'  and  Miles'  mech- 
anism (Barnett,  1968;  Chen  and  Wang,  1983),  i.e., 

G,„(/,d)=  a{f,e)  +  pE{f,8)  (2.47) 

where  a  and  /?  are  linear  and  exponential  growth  coefficients,  respectively.  Contri- 
butions of  the  first  term  (Phillips'  mechanism)  is  small  but  essential  in  their  model 


to  initiate  wave  growth.  The  coefficient  /?  is  proposed  by  Barnett  (1968) 


5s/(^-0.9),    if    ^>0.9 

^  ^  (2.48) 


0,  elsewhere 

Under  this  growth,  wave  energy  Ls  increasing  faster  in  the  high  frequency  compo- 
nents than  the  low  frequency  components.  Therefore,  saturation  condition  needs  to 
be  invoked  to  exclude  the  excessive  energy  by  means  of  wave  breaking.  In  shallow 
water  region,  the  saturation  condition  has  been  suggested  by  Bouws  et  al.(1983) 
to  be  the  product  of  the  deepwater  saturated  spectrum,  such  as  the  one  shown  in 
Eq.(2.34),  and  an  energy  transfer  function 

*<*•'''  -  k'(f,k)c,(f,k)  -  Tnl^UjT)  P-«) 

The  growth  of  a  spectral  component  according  to  the  combined  Phillips'  and 
Miles'  mechanism  can  be  studied  for  the  hypothesized  fetch-limited  and  duration- 
limited  cases.    For  a  steady  state,  one-dimensional,  and  frictionless  situation,  the 


29 

growth  of  a  spectral  frequency  component  in  the  downwind  direction  is  governed 
by  the  following  ordinary  differential  equation 

C,{f)^E{f)  =  a  +  /3E{f) 

By  specifying  E{f)  =  0  at  x  =  0  as  the  boundary  condition,  the  solution  for  spectral 
component  E{f)  is 

For  small  fetch,  E{f)  approximates  {a/Cg)x,  which  corresponds  to  the  Phillips' 
linear  growth  mechajiism.  For  large  fetch,  E{f)  approximates  (a//3)  exp(/3x/Cj), 
which  corresponds  to  the  Miles'  exponential  growth  mechanism.  In  the  real  situa- 
tion the  exponential  growth  is  limited  to  a  saturation  condition  and  an  overshoot 
growth  occurs  between  the  exponential  growth  region  and  the  saturation  region.  A 
typical  growth  curve  for  one  spectral  component  in  the  fetch-limited  case  appears 
m  Fig.  2.4(a). 

For  a  duration-dependent,  one-dimensional,  and  frictionless  situation,  the  growth 
of  a  spectral  component  can  be  determined  from  the  following  ordinary  differential 
equation 

A^(/)  =  a  +  (3E{f) 

By  specifing  E{f)  =  0  at  t  =  0  as  the  initial  condition,  the  general  solution  for  E{f) 

IS 

E{f)  =  ^{expm  -1] 

For  small  duration,  E{f)  approximates  at,  which  corresponds  to  the  Phillips'  mech- 
anism. For  large  duration,  E{f)  approximates  {a//3)  exp{/3t),  which  corresponds  to 
the  Miles'  mechanism.  In  the  real  situation  the  exponential  growth  is  again  limited 
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Figure  2.4:  Schematic  illustration  of  growth  of  a  spectral  component,  £"(/),  at  (a) 
fetch-limited  case  and  at  (b)  duration- limited  case. 
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to  a  sat;iration  condition  and  an  overshoot  growth  occurs  between  the  exponen- 
tial growth  and  the  saturation  periods.  A  typical  growth  curve  for  one  spectral 
component  in  the  diiration-limited  case  appears  in  Fig.  2.4(b). 

The  rate  of  dissipation  in  the  shallow  water,  G^ia{f,6),  is  in  general  thought 
to  come  from,  five  sources:  bottom  friction,  percolation,  whitecaps,  breaking  in 
maintaining  the  equilibrium  spectrum,  ajid  breaking  due  to  the  depth  limitation. 

The  dissipation  rate  attributed  to  bottom  friction  can  be  determined  by  the 
work  done  by  the  bottom  stress,  i.e.. 


Gdi,,i{f,6)  dfdd  =  dndui{f,9)/{p^g) 

where  uj  is  horizontal  velocity  component  at  sea  bed,  r^  is  the  bottom  shear  stress, 
ajid  p^  is  the  water  density.  By  expressing  the  bottom  shear  stress  based  on  the 
empirical  law 

Tb  =  -Pw  Cf  Ub  (/,  0)  \Ub  (/,  6)  I 

where  c/  is  bottom  friction  coefficient,  and  expressing  dui,{f,9)  in  terms  of  dr]{f,0) 
according  to  linear  wave  theory, 

the  rate  Gdi,^i,{f,0)  may  be  approximated  by  the  equation  (Collins,  1972) 

where  (uj)  denotes  the  root  mean  square  value  of  horizontal  wave  orbital  velocity 
component  at  bed  and  can  be  computed  from  the  expression 

=  r  r 'T  "[{■'!  M6f  (2.5,) 

Jo     y-T  a^cosh^  kh  ^        ' 
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The  rate  of  dissipation  of  wave  energy  in  the  shallow  water  due  to  effect  of  a 
porous  sandy  bottom,  Gdis,p{f,  &),  caJi  he  estimated  by  the  work  done  by  the  action 
of  wave- induced  dynamic  pressure  at  the  surface  of  the  bed  to  cause  the  fluid  motion 
in  the  porous  sand  layer. 


Gd.-..p(/,^)  dfdO  =  dp.dw,  l{p^g)  (2.52) 

where  pj,  and  w^,  respectively,  stajid  for  the  dynamic  pressure  and  vertical  water 
partical  velocity  at  the  surface  of  the  bed.  For  porous  sand  layer  of  depth  d  this 
rate  of  dissipation  derived  by  Shemdin  et  al.(l978)  is 


I tanh(^/A^/A2  kd) 

where  Ai  and  A2  are  the  horizontal  and  vertical  percolation  coefficients,  respectively, 
and  they  have  the  dimension  of  (length) /(time). 

The  dissipation  of  wave  energy  caused  by  damping  due  to  kinematic  viscosity, 
z/jt,  and  turbulent  eddy  viscosity,  i>t,  according  to  a  theoretical  analysis  made  by 
Lamb  (1932),  can  be  shown  as 

GdisM^Q)  =  -^[^k  +  Ut)k''E[f,d)  (2.54) 

The  kinematic  viscosity  of  sea  water  i/jt,  which  is  approximately  equal  to  0.018 
cm^/sec  at  32"^,  is  generally  several  orders  of  magnitude  smaller  than  the  turbulent 
eddy  viscosity  Ut.  By  neglecting  the  dissipation  caused  by  the  kinematic  viscous 
effect  in  Eq.(2.54),  and  assuming  the  magnitude  of  turbulent  viscosity  induced  in 
the  water  waves  is  linearly  proportional  to  the  product  of  wave  phase  speed  and 
wave  length,  the  resulting  equation  may  be  expressed  in  the  form 

G4i,Af,&)^-c.fE{f,e)  (2.55) 

where  c„  is  the  constant  of  proportionality.  This  equation  has  been  proposed  by 
Hasselmann  (1974)  to  limit  the  growth  of  wave  energy  on  the  rearface  of  the  spec- 
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According  to  the  JONSWAP  study,  the  dissipation  rate  defined  in  Eq.(2.55) 
is  generally  found  to  be  several  orders  of  magnitude  smaller  than  the  generation 
rate.  However,  actual  dissipation  due  to  whitecaps  may  not  be  small  as  indicated 
by  the  JONSWAP  study  which  utilized  a  hypothesized  equilibrium  spectral  shape, 
i.e.,  JONSWAP  spectrum,  to  define  the  energy  generation  term  that  has  included 
the  influence  of  dissipation  due  to  whitecaps.  Thus,  without  having  a  recisonable 
estimate  of  dissipation  due  to  whitecaps,  the  limit  of  wave  growth  in  the  shallow 
water  model  is  traditionally  given  by  a  saturation  condition  (Chen  and  Wang,  1983): 

A{f)=Bg\27rf)-'^k,h)  (2.56) 

where  B  is  a  coefficient  to  be  of  the  order  of  10~^,  and  #(/:,  h)  is  defined  in  Eq.(2.49). 
The  breaking  of  shallow  water  waves  due  to  depth  limitation  occurs  as  a  conse- 
quence of  strong  wave-bottom  interaction  in  the  nearshore  region.  It  is  this  breaking 
mechanism  that  produces  the  turbulent  flow  and  enhances  the  significant  movement 
of  bed  materials  responsible  for  either  erosion  or  deposition  of  the  bed.  The  estima- 
tion of  energy  dissipation  for  breaking  due  to  depth  limitation  is  important  since  it 
matters  greatly  the  sediment  transport  in  the  cocistal  waters.  However,  an  analyt- 
ical estimation  of  the  dissipation  is  not  found  and  the  determination  of  this  value 
must  rely  on  empirical  formula.  An  empirical  criterion,  which  estimates  the  total 
wave  energy  retained  in  waves  after  the  waves  break  due  to  depth  limitation,  often 
used  in  the  shallow  water  model  is  (Chen  and  Wang,  1983) 

{Etotai,ty''  =  ^  tanh(A:,A)  (2.57) 

where  the  subscript  6  indicates  the  breaking  condition.  The  value  ki,  is  usually 
taken  as  the  wave  number  corresponding  to  modal  period  although  it  is  in  general 
not  clear  what  value  of  ki,  should  be  used  in  the  equation. 
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2.7     Wave  Data  and  Comparisons 

In  1977,  the  Department  of  Coastal  and  Oceanographic  Engineering  at  Univer- 
sity of  Florida  (UF)  initiated  a  shallow  water  wave  data  collection  program  around 
the  State  of  Florida.  By  1985,  twelve  wave  stations  have  been  established  (See 
Fig.  2.5).  In  the  offshore  region,  there  are  two  NOAA  maintained  wave  buoys  with 
ID  No.  41002  and  41006  as  shown. 

Of  the  twelve  shallow  water  stations,  data  from  Kings  Bay,  Marineland  and 
Cape  Caneveral  are  eventually  being  used  here  for  comparison  purposes.  At  Kings 
Bay,  two  PUV  imderwater  packages  as  majiufactured  by  Sea  Data  Inc.  (Model 
635-12)  were  deployed  at  15  m  and  18  m  water  depth.  These  packages  are  self- 
contained  with  the  ability  to  measure  two  axes  water  particle  velocities  U  and  V, 
and  the  pressure  information,  P.  Directional  wave  information  is  derived  form 
the  P-U-V  data  based  on  the  method  proposed  by  Longuet-Higgins  et  al.(l963). 
Detailed  data  handling  ajid  ajialysis  procedures  are  given  in  the  User's  Guide  by 
UF  (1984).  The  gages  at  Marineland  aoid  Cape  Caneveral  are  shore-connected 
imderwater  pressure  transducers  at  10  m  and  8  m  water  depths  respectively.  Only 
wave  height  information  can  be  retrieved  from  them. 

The  deepwater  buoy  stations  are  located  at  29.3''N,  77.3°W  (#41006)  and  at 
32.3''N,  75.3° W  (#41002).  The  water  depths  of  these  two  stations  are  970  m  and 
3200  m,  respectively.  These  buoys  also  collect  5  m-level  winds  every  hour,  baised  on 
aji  8.5  minute  record. 

One  of  the  intentions  of  the  field  data  measurement  program  is  to  provide  wave 
climate  information  along  the  south  eeistern  seaboard  of  the  United  States.  To 
accomplish  this  goal,  a  numerical  wave  prediction  model  is  required  to  fill  in  the 
data  gaps  and  to  extend  data  extrapolation. 

The  WIS  model  was  first  tested  against  the  measured  deepwater  wave  data 
at  NOAA  buoy  #41006.    A  segment  of  the  data  covering  the  complete  month  of 
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Figure  2.5:  Wave  gage  stations. 
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Figure  2.6:  (a)  Weather  map  showing  northeajst  winds  generated  by  a  high  pressure 
system  (1/12/84). 

January,  1984  was  used  for  testing  purposes  for  two  reasons:  (i)  During  this  month, 
directional  wave  data  at  Kings  Bay  was  available  for  the  entire  month  and  (ii) 
January  is  usually  the  month  of  high  waves  dominated  by  large  scale  high  pressure 
systems.  Figure  2.6(a)  and  2.6(b)  show  two  examples  of  typical  large  scale  high 
pressure  systems  occxirring  in  January.  In  these  examples,  the  northeast  winds 
developed  from  the  high  pressure  systems  will  generate  waves  that  propagate  to  the 
southeastern  coast  of  the  United  States. 

The  grid  system  covers  an  orthogonal-spherical  area  with  the  east  and  the  west 
boundaries  coinciding  with  70.0''W  and  82.5''W,  respectively,  and  with  the  north 
and  south  boundaries  coinciding  with  42.5''N  and  22.5''N,  respectively.  The  grids  are 
all  parallel  to  the  longitudes  and  latitudes  such  that  each  grid  is  2.5  degrees  apart 
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Figure  2.6:  (b)  Weather  map  showing  northeast  winds  generated  by  a  high  pressure 
system  (1/21/84). 
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Figure  2.7:  Deepwater  grid  coordinate  system. 


from  the  next  and  a  total  of  6x9  grid  points  are  found  in  the  system  (Fig.  2.7). 
The  wind  information  served  as  the  input  for  the  model  was  obtained  by  analyzing 
surface  pressures  from  daily  weather  maps.  The  theories  of  the  geostrophic  winds 
and  plajiet  boundary  layers  were  used  in  the  analysis  (see  Section  3.4  -  wind  model). 
The  10  m-level  winds  were  generated  at  the  grid  system  for  forty  days  starting  on 
December  21,  1983.  The  computed  winds  were  then  blended  with  the  ones  measured 
by  the  NOAA  buoys. 

The  radiative  boundary  condition  is  utilized  in  the  wave  model.  However,  two 
modifications  were  made  to  the  boundary  condition.  The  first  one  was  made  to 
the  boundary  where  the  grid  system  is  connected  with  land.    It  simply  regulates 
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that  no  wave  energy  is  transported  to  sea  from  land.  The  second  one  was  made 
to  the  boundary  where  the  grid  system  is  connected  with  the  sea.  It  regulates 
that  the  amount  of  wave  energy  advected  from  a  grid  point  at  the  boundary  to 
the  grid  points  Inside  the  boundary  is  the  same  as  the  wave  energy  transported 
from  the  area  outside  the  boundary  to  the  grid  point  at  the  boundary.  Therefore, 
a  quasi-stationairy  wave  energy  transport  process  is  assumed  at  the  seaboard  grid 
boundaries. 

The  wave  model  was  run  at  time  increment  equal  to  two  hours  for  a  total  of 
forty-one  days  from  December  21,  1983  to  January  31,  1984.  The  first  ten  days  were 
for  initialization  and  the  results  were  not  used  for  comparison.  Figxire  2.8  shows  the 
comparison  of  the  time  series  of  the  wave  energy  spectral  evolution  process  of  the 
entire  month.  The  correlations  of  the  mezisured  and  computed  modal  period  and 
significant  wave  height  are  shown  in  Fig.  2.9.  The  significant  wave  height,  which  is 
commonly  denoted  as  Hijz  or  H^,  is  defined  as  the  average  height  of  the  one-third 
highest  waves.  When  the  heights  of  individual  waves  in  a  deepwater  wave  record 
are  ranked  from  the  highest  to  lowest,  the  frequency  of  occurrence  of  waves  above 
any  given  value  is  given  to  a  close  approximation  by  the  cumulative  form  of  the 
Rayleigh  distribution.  Therefore,  the  significant  wave  height  can  be  determined 
approximately  from  the  following  equation: 

H,  =  4.0  (Et^taiV^' 

The  correlation  coefficient,  6,  given  in  Fig.  2.9  is  defined  as  follows: 

6  =  t,li^MX^),-{XMXrn)i]/{[h{Xc)i-mi]'hi^rn)i-lX^^  (2.58) 

1=1  1=1  «=1 

where  Xc  is  a  computed  X,  Xm  is  a  measured  X,  N  is  the  sample  size,  and  the 
overhead  bar  indicates  the  operation  of  arithmatic  meaji. 

The  measured  and  hindcztsted  energy  scale  parameter,  f^H,/g,  are  plotted  in 
Fig.  2.10  for  waves  with  H,  >3  m,  3  m>  H,  >  2  m,  and  H,  <2  m. 
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Figure  2.8:   Comparison  of  NOAA  and  WIS  wave  spectra  for  January  1984  at  the 
Buoy  #41006  location. 
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Figure  2.9:    Comparisons  of  NOAA  and  WIS  modal  periods  and  significant  wave 
heights  at  the  Buoy  #41006  location. 
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Figure  2.10:  Comparison  of  NOAA  and  WIS  wave  energy  scale  parameters  at  the 
Buoy  #41006  location. 
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Clearly,  the  model  performs  better  for  hindcasting  wave  height  (  with  a  cor- 
relation of  0.5)  than  for  wave  frequency  (with  a  negative  correlation  coefficient). 
The  model  performance  also  improves  progressively  for  higher  wave  environment 
(Fig.  2.9)  and  this  is  particularly  the  case  if  the  wind  wave  dominates.  The  model 
faxes  the  worst  on  hindcasting  the  swell  component.  As  can  be  seen  in  Fig.  2.8, 
both  growth  and  decay  of  low  frequency  components  hindcasted  by  the  model  are 
clearly  inadequate. 


CHAPTER  3 
A  UNIFIED  COASTAL  WAVE  PREDICTION  MODEL 


3.1     Introduction 

It  is  shown  in  the  last  chapter  that  the  parametric  model  and  the  discrete  model 
are  two  alternative  methods  to  hindcast  waves  in  deep  water.  Growth  rates  of  wave 
height  predicted  by  the  parametric  model  of  Ha^selmaim  et  al.(1976)  and  the  WIS 
model  were  consistent  in  the  case  which  waves  are  growing  either  along  a  fetch  or 
through  time  (Resio,  1981;  Resio  and  Tracy,  1983).  However,  the  comparison  of  the 
hindcast  based  upon  either  model  to  the  measured  wave  heights  from  buoy  in  the 
open  oceaji  often  shows  discernible  deviations  (Gunther  et  al.,  1979;  Corson  and 
Resio,  1981).  The  comparison  of  the  modal  periods  computed  from  either  model 
to  the  buoy  data  also  shows  certain  degree  of  discrepancy  (Gunther  et  al.,  1979). 
Thus,  it  seems  necessary  to  modify  the  existing  models  in  order  to  improve  the 
accuracy  of  the  hindcasted  waves  compared  to  the  measured  data. 

A  new  unified  wave  model,  named  here  as  the  Unified  Coastal  Wave  Predic- 
tion (UCWP)  model,  was  developed  in  the  present  study.  The  model  is  of  discrete 
type  and  is  comprised  of  deepwater  and  shallow  water  submodels.  The  deepwater 
submodel  is  a  modified  WIS  model.  The  modifications  ase  mainly  in  the  source 
fimctions  by  incorporating  some  recent  development  in  air-sea  interaction.  The 
model  is  calibrated  by  utilizing  the  field  wind  and  wave  data  collected  by  NOAA 
buoys  ajid  it  is  suitable  for  hindcasting  ocean  waves  generated  by  either  high  pres- 
sure or  low  pressure  weather  system.  The  shallow  water  submodel  is  baised  upon 
the  numerical  algorithm  developed  by  Chen  and  Wang  (1983)  and  it  is  suitable  for 
hindcasting  coastal  water  waves. 
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In  the  following  sections  of  this  chapter,  both  the  UCWP  deepwater  and  shallow 
water  submodels  will  be  introduced.  The  deepwater  submodel  is  to  be  described  in 
Section  3.2  and  the  shallow  water  submodel  is  in  Section  3.3.  A  wind  model  which 
prepares  the  wind  data  for  the  generation  of  water  surface  waves  in  the  UCWP 
model  Ls  described  in  Section  3.4. 

3.2     UCWP  Deepwater  Submodel 

The  deepwater  submodel  is  a  spectral  model  of  the  discrete  type.  This  model 
is  more  flexible  in  selecting  source  functions  than  a  pareimetric  model  as  it  is  not 
always  aoi  easy  task  to  introduce  new  source  functions  in  the  JONSWAP  parameter 
space. 

The  governing  equations  for  the  UCWP  deepwater  model  are  the  same  as  the 
WIS  model  described  earlier,  i.e.,  Eqs.(2.27)  ajid  (2.28).  The  advection  term  in 
model  is  computed  in  the  same  manner  as  in  the  WIS  model.  The  source  terms 
including  both  generation  and  dissipation  in  the  UCWP  deepwater  model  are  mod- 
ified from.  WIS  model. 
3.2.1     Generation  terms 

Recent  studies  made  by  Mitsuyasu  and  Honda  (1982),  and  Mitsuyasu  (1985) 
indicate  that  the  mechanism  corresponding  to  wave  growth  due  to  wind  input  shall 
be  considered  separately  for  the  prevailing  gravity  waves  and  the  locally  generated 
wind  waves.  Therefore,  the  utilization  of  Eq.(2.23)  alone  to  estimate  wave  growth 
does  not  seem  to  be  a  proper  approach  in  the  parametric  model.  On  the  other  hand, 
the  WIS  model  which  utilizes  two  separate  mechanisms  expressed  by  Eqs.(2.29)  and 
(2.31)  to  estimate,  respectively,  the  growth  of  prevailing  waves  cind  new  wind  waves 
seems  to  be  more  consistent  with  the  nature  of  wave  growth. 

In  the  WIS  model,  the  wave  generation  comprises  of  three  terms  -  growth  of 
prevailing  waves,  generation  of  new  wind  waves  and  nonlinear  trainsfer  of  energy 
among  wave  components.  The  growth  terms  are  computed  at  time  level  k  +  1  and  the 
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transfer  term  is  computed  at  time  level  k.  The  first  ajid  the  last  terms  -  Eqs.(2.29) 
and  (2.43),  respectively  -  are  of  a  sim^ilar  nature,  responsible  for  the  evolution  of 
existing  wave  components.  Their  combined  effect  is  shown  in  Fig.  3.1.  Here,  the 
nonlinear  energy  transfer  is  seen  to  be  the  dominant  mechanism  and  has  an  abrupt 
cutoff  at  /  =  /m.  The  growth  of  the  down-slope  face  of  the  spectnma  is  shown  to 
be  a  minor  process.  At  peak  spectral  frequency,  /  =  fm,  approximately  80  percent 
of  the  energy  growth  is  being  transferred  to  lower  frequencies.  By  comparison  with 
field  data  as  discussed  in  the  previous  section,  the  predicted  spectral  shape  is  found 
to  be  too  sharp  peaked  due  to  the  abrupt  cutoff  aind  the  low  frequency  components 
grow  too  fcist  because  of  the  proposed  nonlinear  transfer  process.  A  new  growth 
terni  is  proposed. 

Since  we  have  shown  that  the  wave-wave  interaction  source  term  behaves  like 
an  apparent  wind  source  term,  it  can  be  lumped  into  the  wind  source  term  for 
the  growth  of  the  prevailing  waves.  Based  upon  the  common  practice  that  this 
source  term  should  be  related  to  a  stable  spectral  shape  to  maintain  similarity, 
we  confine  our  search  to  established  spectral  functions  by  testing  them  against  field 
data  collected  during  steady  winds  when  the  prevailing  wave  growth  is  the  dominant 
mechanism.  After  various  trails,  a  new  source  function  is  proposed  here: 


Gin,x{fJ)   = 


c,a'fEjM{f)^cos\9-d),      if|5-5|<^   and   |^,-^|<^ 

oTT  Z  2 


0,  elsewhere 

(3.1) 

where  c,  is  a  non-dimensional  coefficient,  regarded  by  some  investigators  to  be  a 

function  of  density  ratio  of  air  to  water,  and  EjM{f)  is  the  one-dimensional  modified 

JONSWAP  spectrum  proposed  recently  by  Donelan  et  al.(1985): 

EjM)  =  j^^^^^exp[-{Ufy]^  'i'o^fj^  (3.2) 

with  all  its  free  parameters  a,  /^,  7,  and  a^t  defined  the  same  as  in  the  JONSWAP 
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Figure  3.1:  Combined  sources  of  nonlinear  term  and  generation  term  in  WIS  model. 
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spectnim.  Substituting  Eq.(3.2)  into  Eq.(3.1)  and  absorbing  c,  into  a,  we  have 

IJ-fr, 


GinAfJ)={ 


if  \0-6\<-  and  |^^  -0\<- 
0,  elsewhere 


(3.3) 


The  peak  frequency  j^  is  determined  from  the  following  equation: 

!):'   =  {UL)-'"'  +  ^^(^10COs(^-^J)4/3^,^-3/7  (34) 

which  is  similar  to  Eq.(2.39)  with  the  exceptions  that  the  wind  component  parallel 
to  the  mean  wave  propagation  is  used  here  instead  of  the  absolute  strength  and 
that  the  spectral  peak  frequency  /^  given  in  Eq.(2.39)  is  now  replaced  by  a  newly 
defined  peak  energy  transfer  frequency,  fma-  This  peak  energy  transfer  frequency  is 
defined  as 

^(/^,)=  max  {£:{/)  >^,(/)} 

where  Edf)  is  a  threshold  spectrum  below  which  no  energy  transfer  will  occur  from 
high  frequency  components  to  low  frequency  components.  To  use  fma  instead  of  fm 
in  Eq.(3.4)  physically  limits  the  energy  transfer  to  those  components  exceeding  the 
threshold  condition. 

At  present,  there  is  no  established  guideline  to  determine  the  threshold  spec- 
trum, Ec{f).  However,  based  upon  physical  reasonings,  we  observe  here  that  £'c(/) 
should,  at  least,  satisfy  the  following  criteria:  (i)  The  threshold  spectrum  should 
fall  below  the  saturation  (equalibrium)  spectrum  eis  energy  transfer  is  expected  to 
take  place  before  the  spectral  energy  reaches  saturation,  (ii)  The  threshold  spec- 
trum should  be  a  monotonically  increasing  function  with  decreasing  frequency  and 
should  approach  equalibrium  spectrum  at  low  frequencies  where  energy  transfers 
become  more  and  more  difficult. 
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To  propose  a  suitable  Edf),  two  years  measured  data  -  1984,  1985  -  from 
NOAA  buoys  #41002  and  #41006  were  examined.  Since  we  are  dealing  with  the 
mechanism  on  the  growth  of  prevailing  waves,  only  those  data  corresponding  to 
increasing  in  both  wave  energy  ajid  modal  period  were  used.  Figure  3.2  shows  the 
plots  of  E{fm)  versus  fm  from  the  data  of  buoys  #41002  and  #41006.  If  the  mean 
values  of  the  data  points  are  connected,  they  form  a  skewed  bell-shaped  curve  with 
a  peak  value,  fp,  occurring  at  frequency  between  0.09  and  0.11  Hz.  This  indicates 
that  energy  transfer  experiences  strong  resistance  when  fm  becomes  less  than  this 
peak  value. 

Based  upon  the  above  observations,  we  propose  here  the  following  functional 
form  for  Edf): 

Ecif)  =  <=ms-0fj-,e-'^'exv[-{j-r]  (3.5) 

which  is  the  P-M  equilibrium  spectrtmi  modified  by  a  multiplier.  Here  Cms  is  an 
empirical  scale  factor  which  is  chosen  as  equal  to  0.183  by  matching  Eq.(3.5)  with  the 
data  mean  at  designated  /p=0.1  Hz.  Both  Eq.(3.5)  and  the  original  P-M  equilibrium 
spectrxim  are  shown  in  Fig.  3.2,  with  the  data.  As  can  be  seen,  Ec{f)  has  the  basic 
characteristics  as  discussed  with  the  most  active  energy  transfer  process  taking  place 
in  the  energy  containing  range. 

It  is  noted  here  that  plots  of  spectral  peak  energy  versus  spectral  peak  frequency 
in  the  wave  growth  cases  in  Fig.  3.2  show  that  the  data  are  scattered  in  the  equi- 
librium range  indicated  by  the  P-M  spectrum.  Although  the  data  scattering  could 
be  exaggerated  by  the  presence  of  swells,  which  generally  appear  as  long  waves, 
there  is  no  indication  that  a  simple  curve,  such  as  the  P-M  equilibrium  spectrum, 
can  represent  the  relationship.  The  data  scattering  suggests  active  nonlinear  wave- 
wave  interactions.  That  is,  energy  transfer  from  higher  frequency  components  to 
lower  frequency  ones  is  a  nonlinear  process  and  it  is  dependent  upon  multi-variables 
including  wind  speed,  wave  phase  speed,  peak  energy  frequency,  etc. 
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Figure  3.2:  Plots  of  E{fJi  versus  /„ 
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The  nonlinear  transfer  coefficient,  a„i,  in  Eq.(3.4)  has  been  suggested  by  Has- 
selmann  (1976)  to  assume  a  constant  value  of  0.00138  under  the  condition  of  steady 
wind  within  a  duration  At.  However,  the  values  of  a„/  evaluated  from  the  present 
NOAA  buoy  data  appear  to  be  varying  instead  of  constant.  Figure  3.3  shows  a^ 
values  from  buoys  #41002  and  #41006  plotted  against  the  non-dimensional  param- 
eter fo/{f^^At).  It  is  clearly  seen  there  exists  a  lower  bound  of  the  plotted  data 
which  are  approximated  by  the  dashlines  in  the  figure.  Above  the  dashlines,  most 
data  points  are  clustered  in  a  narrow  strip  shown  in  the  thickline  box.  Inside  the 
box,  the  separation  of  data  points  into  two  large  clusters  is  caused  by  the  shortage 
in  having  the  one-hour  time  interval,  At,  in  computing  ani  in  Fig.  3.3.  Outside 
the  box,  the  data  shows  considerable  scatter,  which  is  due  to  the  contamination  of 
swells  moving  into  the  region  (data  with  dual  spectral  peaJcs).  If  a  curve  is  fitted 
to  the  upper  bound  of  the  clustered  data  and  using  the  condition  that  the  value  of 
a„i  approaches  the  Hasseimann's  constant  of  0.00138  when  fo/{f^^At)  —>■  0  (steady 
wind  over  long  period),  we  arrive  at 

aru  =  26(0.0195  +  i-^/^'  f  (3.6) 

This  is  drawn  in  Fig.  3.3  as  a  solidline  curve  which  approaches  asymptotically  to  the 
dashline,  the  lower  bound  of  the  plotted  data,  as  fo/{f^^At)  approaches  infinity. 
Equation  (3.6)  has  been  used  in  the  new  model  as  a  unique  parametric  formula  in 
accounting  the  migration  rate  of  spectral  peak  in  the  frequency  domain. 

The  same  data  set  excluding  those  exhibit  the  swell  property  is  then  employed 
to  determine  the  parametric  formulas  for  the  parameters  a,  7,  a^,  and  at,.  This  is 
accomplished  by  minimizing  the  square  errors  over  the  ensemble  of  measured  data, 
AE{f)/At,  and  the  computed  ones  from  the  proposed  Eq.(3.3).  The  values  of  a 
are  plotted  against  fmi/ fo  in  Fig.  3.4.  The  data  are  seen  more  scattering  when 
fmt/fo  is  less  than  0.5,  which  indicates  that  the  generation  of  local  wind  waves 
becomes  important  in  this  situation  to  effect  the  nonlinear  transfer  process.   Since 
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the  values  of  AE{f)/At  from  the  data  are  contaminated  by  the  advected  wave 

energy  into  the  region,  the  best-fit  curve  of  the  lower  bound  values  is  believed  to  be 
more  representative  of  the  growth  rate.  This  curve  can  be  approximated  by 

a  =  0.0122  (^)'/^  (3.7) 

/o 

or  equivalent  ly, 

a  =  0.0122  {^y/^ 

where  Cms  =  gl2-K f^t  is  the  wave  celerity  at  fma-  The  ratio  fmal fa  is  n«imed  the 
resonant  parameter,  which  has  been  found  important  to  estimate  the  nonlinear 
transfers  of  energy  in  previously  developed  models,  for  instance,  it  appeared  in 
Eqs.(2.36)  and  (2.37),  respectively,  in  the  parametric  and  WIS  models.  Also  shown 
in  Fig.  3.4  axe  the  empirical  equations  proposed  for  a  by  Hasselmann  et  al.  (1976) 
cind  Donelan  et  al.(l985).  The  equation  proposed  by  Hasselmann  et  al.  (1976)  is 
obtained  by  defining  a  constant  wind  speed  during  Ai  in  Eq.(2.36),  which  is  derived 
in  the  parametric  model  by  allowing  the  variations  of  only  a  and  fm- 

a  =  0.0091(^)2/^ 
/o 

The  empirical  equation  for  a  proposed  by  Donelan  et  al.(1983)  is 

a  =  0.006(^)°" 

/o 

which  is  obtained  based  upon  data  collected  in  the  Great  Lakes.  It  should  be  noted 
here  that  when  the  above  two  equations  were  compared  with  Eq.(3.7)  in  Fig.  3.4 
they  were  multiplied  by  a  factor  (4.4)^/^  which  is  the  constant  absorbed  in  a  in 
Eq.(3.7).  Both  empirical  equations  proposed  by  Hasselmann  et  al.  (1976)  and 
Donelan  et  al.  (1985),  as  shown  in  Fig.  3.4,  appeared  to  be  less  representative  for 
the  lower  boimd  values.  The  equation  proposed  by  Hasselmann  et  al.  (1976)  is 
more  representative  for  the  mean  values.  In  the  present  study,  Eq.(3.7)  is  adopted. 
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Figure  3.5:  Plot  of  -7  computed  from  Eq.(3.8)  versus  7  from  the  least  square  method. 


The  7  values  solved  from  the  least  square  method  can  be  approximated  by  the 
following  equation  (see  Fig.  3.5): 


7  =  exp[  2 


/m,|  ^2  5  (  Am^3 


(3.8) 


fo"        Vo 

The  above  equation  gives  a  minimum  7  «  2.38  when  fms/  fo  =1-5.  Although  7  is 
not  bounded  for  small  fm»l  fo,  the  spectral  peak  value  of  G,„,i,  which  is  proportional 
to  a^7  ^  exp(2/„l,//o)(/„»,//o)^/^  is  bounded  for  small  fm,l  fo- 

Values  <7a  and  Of,  obtained  from  the  field  data  are  scattered  (see  Fig.  3.6).  The 
arithmetic  mean  values  of  Oa  and  oi,  are  0.12  and  0.08,  respectively,  and  they  are 
in  the  same  order  of  magnitude  as  given  by  the  JONSWAP  spectrum.   The  values 
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of  (Ta  and  o-j  control  the  spectral  width  of  the  forward  face  and  down-slope  side  of 
the  spectrum,  respectively.  Since  at  any  instant,  the  energy  contained  In  the  wave 
system,  which  is  manifested  by  the  magnitude  and  shape  of  the  energy  spectrum, 
is  the  cumulative  difference  between  the  gain  and  the  loss.  The  measured  a^  and  C}, 
could  also  be  viewed  as  parameters  characterizing  the  net  energy  levels  on  each  side 
of  the  spectrum.  Since  on  the  down-slope  side  (high  frequency  components),  the 
loss  of  energy  due  to  dissipation  and  transfer  is  expected  to  be  very  active  during 
the  growth  period,  the  measured  ai,  values  evidently  axe  too  small  to  represent  the 
gross  gain.  Based  upon  our  estimate,  a  value  of  approximately  four  times  as  large 
should  be  used,  i.e.,  Of,  should  be  adjusted  to  0.32.  On  the  forward  side,  on  the 
other  hand,  the  loss  due  to  dissipation  is  minimal.  Therefore,  the  arithmetic  mean 
value  of  measured  Oa  should  represent  good  estimate  of  energy  gain  due  to  energy 
trajisfer  and  no  adjustment  is  made. 

The  source  function  G,>,_i  defined  in  Eq.(3.3)  is  also  plotted  in  Fig.  3.1  for 
comparison.  The  energy  growth  on  the  down-slope  face  of  the  spectrum  is  seen  to 
be  much  stronger  thein  that  in  the  WIS  model,  and  the  unnatural  abrupt  cutoff  is 
smoothened  out. 

In  the  calibrations  of  parameters  a  and  7,  the  past  practice  is  to  assume  wind 
direction  has  to  be  parallel  to  the  mean  wave  direction.  However,  in  the  real  situa- 
tion, winds  and  waves  may  have  different  directions,  and  the  wind-wave  resonance 
can  be  greatly  reduced  when  the  wind  direction  deviates  significantly  from  the 
wave  propagation  direction.  Therefore,  this  difference  between  the  wind  and  wave 
directions  needs  to  be  properly  incorporated  into  the  resonance  condition.  This  is 
accomplished  in  the  model  by  redefining  /o  =  g/2irUiocos(d  —  0^)  in  Eqs.(3.7)  and 
(3.8),  which  are  used  to  estimate  a  and  7,  respectively. 

The  rate  of  growth  of  newly  generated  waves  in  WIS  model  is  treated  in  such  a 
way  that  the  total  energy  growth  during  time  step  At  is  independent  of  the  length 
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of  At.  This  requires  either  the  wind  to  be  sufficiently  strong  or  the  duration  to 
be  sufficiently  long,  so  that  in  every  time  step  the  growth  of  local  waves  reaches 
its  upper  limit.  The  WIS  model  was  found  to  have  a  tendency  of  overpredicting 
wave  growth  under  mild  wind  condition  (see  Fig.  2.8)  and  this  source  term  might 
be  partially  responsible. 

We  observe  here  that  the  growth  of  local  waves  due  to  winds  at  the  current  time 
step  can  take  place  in  one  of  the  following  three  forms: 

a.  If  the  shift  in  wind  direction  from  the  previous  time  step  is  sufficiently  large, 
say  over  a  certain  threshold  value,  0c,  new  fresh  waves  will  be  generated  in  the 
direction  of  new  wind. 

b.  If  the  shift  in  wind  direction  from  the  previous  time  step  is  small  and  the 
wind  direction  is  different  from  the  prevailing  wave  direction,  wind  waves  generated 
from  the  previous  time  step  will  continue  to  grow  in  the  direction  of  the  wind. 

c.  If  the  shift  in  wind  direction  from  the  previous  time  step  is  small  and  the 
wind  direction  is  the  same  as  the  prevailing  wave  direction,  this  source  term  will 
vanish  and  growth  will  be  limited  to  prevailing  waves  only. 

Accordingly,  the  source  terms  for  the  generation  of  new  wind  waves,  G,n,2)  are 
proposed  as  follows: 

Case  a.  A^^;  >  9^  In  this  case,  it  is  generally  recognized  that  two  wave  generation 
mechanisms  are  at  work,  the  Phillips'  (1957)  resonance  mechanism  and  the  Miles' 
(1957,  1959a,b,  1962)  instability  mechanism.  It  is  also  recognized  that  for  practical 
application,  the  resonance  mechanism  which  is  only  important  at  the  very  initial 
stage  can  be  neglected.  Therefore,  it  is  common  to  assume  that  the  instability 
mechanism  dominates,  or  G,„,2  oc  E{f,9).  We  further  observe  that  the  growth  at 
a  location  will  eventually  reach  duration  limited  condition,  the  following  functional 
form  should  be  satisfied: 

GinM^O)  =  b{E^{f,d)  -  E{f,9)]  (3.9) 
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where  6  is  a  growth  coefficient  and  Eoo  is  the  duration  unlimited  spectrum,  and  is, 
in  general,  fetch  dependent.  For  the  present  model  of  deepwater  application  b  is 
chosen  as  equal  to  0.1s/(l7io/C)'*  with  s  the  density  ratio  of  air  to  water  and  C  the 
wave  celerity  at  frequency  /,  and  ^oo  is  the  P-M  spectrum  with  the  scale  parameter, 
OLp,  equal  to  0.015,  i.e.. 


Gm.2(/,5)   = 


0.1sf{^y[E^{f)-cos'{9  -  e^)  -  E{f,9)  ],     if  1^  -  ^^1  <  7r/2 

0,  elsewhere 

(3.10) 


The  above  equation  is  applied  in  the  present  model  to  initiate  wave  generation  at 

high  frequency  (/o  >  0.2  Hz)  and  whenever  the  shift  in  wind  direction  exceeds  d^. 

At  present,  there  is  no  known  criterion  on  selecting  dc-   For  convenience,  we  have 

chosen 

^c  =  (— rad/hr)  •  Af  (3.11) 

Is 

Case  b.  A^^,  <  0c  In  this  case,  the  growth  mechanism  should  be  similar  to  that  of 
the  prevailing  wave  components  as  discussed  earlier  with  the  exception  that  now  the 
growth  should  be  centered  around  the  wind  direction,  0^,  instead  of  the  prevailing 
wind  direction,  0.  Thus,  the  source  function  is  aissumed  to  have  similar  form  as 
Eq.(3.3)  with  0  being  substituted  by  0^,  i.e.. 


'm,2 


{f,0) 


(4  9^        r     Jmo.^y    8 


exp[-(^)']  —  cos^(5-^^),       if  \0-0^\<-kI2 


(27r)VVmo      "^     ^   /   '  ^  37r 

0,  elsewhere 


(3.12) 
The  spectral  peak  frequency  f^o  and  scale  parameter  a  are  determined  by  the  Scime 
equations  used  in  the  prevailing  wave  case,  i.e.,  Eqs.(3.4)  and  (3.7),  or 


and 

/o 


a,,  =  0.0122  (=^)^/^  (3.14) 
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with  /o  =  g/2irUioCOs{A6^)  being  defined  in  Eq.(3.14)  to  account  for  the  influence 

of  Ad^  to  the  resonajice  condition. 

Case  c.  9^  =  9  In  this  case,  the  wind  direction  is  the  same  as  the  prevailing  wave 
direction,  we  have  Gm,2(/i^)  =  0- 

Particular  interest  should  be  noted  to  the  source  function  given  in  Eq.(3.10).  If 
this  is  the  only  source  term  for  spectral  evolution  at  the  duration-limited  case,  the 
transport  equation  can  be  solved  ajialytically  to  yield  the  spectrum 

E{f,9)  ^  E^{f){l  -  exp[-0.1sf{^)U]}-cos'{9  -  9^)  (3.15) 

which  vanishes  at  i  =  0,  and  is  equal  to  E^{f)  [2 / n)  cos'^ [9  —  9yj)  at  f  =  oo.  A 
simulation  of  the  growth  of  wave  spectrum  based  on  Eq.(3.15)  under  30-knot  winds 
is  shown  in  Fig.  3.7.  The  result  is  compared  to  Barnett's  (1968)  model,  also  shown  in 
Fig.  3.7,  and  the  comparison  shows  a  good  agreement  between  the  two  simulations 
although  the  spectra  predicted  based  on  Eq.(3.15)  have  wider  and  flatter  shape  than 
these  predicted  from  Barnett's  model. 

If  Eoo{f)  is  eliminated  from  Eq.(3.10)  by  using  Eq.(3.15),  the  source  function 
Ginaiff^)  ^^^  ^®  rewritten  as 

This  equation  indicates  that  the  growth  rate  is  equal  to  the  product  of  a  spectrum 
ajid  a  forcing  function  which  is  dependent  on  wind  strength,  wave  frequency,  and 
time.  The  relation  described  in  Eq.(3.16)  is  diff"erent  from  Barnett's  formula  de- 
scribed in  Eq.(2.48)  which  states  that  the  growth  rate  is  equal  to  the  product  of 
a  spectrum  and  a  forcing  function  which  is  dependent  only  on  wind  strength  and 
wave  frequency.  Since  growth  rate  proposed  by  Barnett  (1968)  can  not  predict  a 
steady  state  spectrum  and  it  has  an  artificial  cutoff"  of  energy  growth  as  the  ratio 
Uio/C  is  less  than  0.9,  it  is  more  appropriate  to  use  Eq.(3.10)  rather  than  Barnett's 
formula  in  predicting  the  local  generation  of  wind  waves. 
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Figure  3.7:  Simulations  of  spectral  growth  at  duration-limited  case;  U\q  =  30  knots. 
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3.2.2     Dissipation  Terms 

The  dissipation  of  wave  energy  caused  by  turbulence  heis  been  studied  by  a  num- 
ber of  investigators.  In  deep  water,  the  functional  form  suggested  by  Hasselmann 
(1974)  is  represented  by 

G,i.{f,e)  =  -c.fE[f,9)  (3.17) 

where  c„  is  an  dissipation  coefficient  and  its  value  is  usually  obtained  by  utilizing 
field  wave  data  obtained  during  successively  decreasing  winds  with  steady  direction. 
The  value  selected  in  the  present  model  is  s/6.  Substituting  this  value  into  the  above 
equation  and  integrating  at  the  duration-limited  case,  we  obtain 

E\l\f,9)  =  El,{f,e)  e-^f^'l^  (3.18) 

Figure  3.8  offers  a  comparison  of  wave  decays  between  computed  values  and  field 
data.  Owing  to  this  calibration  practice,  the  above  dissipation  function  is  found  to 
be  more  suitable  for  laxge  weather  systems  such  as  related  to  high  pressure  cells,  in 
which  the  wind  direction  is  relatively  stable. 

For  low  pressure  systems,  such  as  hurricane  winds,  the  wind  direction  could  be 
completely  opposite  on  the  forweird  cind  the  trailing  phases  of  the  eye  while  the 
wind  strengths  are  usually  on  the  same  order  of  magnitude.  With  such  a  weather 
system  the  process  of  wave  energy  dissipation  is  much  more  vigorous  as  wind  in  the 
trailing  phase  often  blows  against  the  prevailing  waves  and,  thus,  promotes  wave 
breakings.  For  low  pressure  systems,  the  dissipation  function  may  be  determined 
based  on  that  proposed  by  Lamb  (1932): 

Giu{f,e)  =  -c,fE{f,0)  -  -M_E{f,e)  (3.19) 

Pu,C  |dr/| 

where  p  is  the  normal  pressure  disturbance  and  rj  is  the  water  surface  fluctuation.  By 
assuming  dp  is  of  the  same  order  of  magnitude  as  the  downwind  flux  of  the  horizontal 
momentum  of  the  wind,  d{paU^o)  (Phillips,  1967;  Longuet-Higgins,  Cartwright,  and 
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Figure  3.8:  Comparison  of  spectral  decays  between  measured  and  computed  values. 
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Smith,  1963),  Eq.{3.19)  can  be  approximated  by  the  following  expression: 

G,is{f,9)  =  -c,fE{f,d)-Qm^E{f,d)  (3.20) 

which  is  a  modification  of  that  in  high  pressure  systems  by  adding  an  additional  term 
—O.Q2[sUlQJCH,\E{f,9).  K  this  is  the  only  soiirce  term  in  the  transport  equation, 
the  wave  energy  evolution  at  the  duration-limited  condition  can  be  found  as 

where  A  =  0.02sU^q/ CH,.  The  charateristics  of  e~^^'  are  such  that  it  varies  from  0 
to  1  depending  upon  the  maturity  of  the  sea  and  that  smaller  waves  dissipate  more 
energy  against  wind  than  that  of  larger  waves. 

3.3     UCWP  Shallow  Water  Submodel 

The  shallow  water  submodel  is  based  upon  the  numerical  algorithm  developed 
by  Chen  and  Wang  (1983)  which  is  a  nonstationary  spectral  transformation  model  of 
discrete  type.  The  governing  equations  utilized  in  the  model  are  given  in  Eqs.(2.45) 
and  (2.46).  In  the  present  application,  a  number  of  source  terms  are  modified. 
3.3.1     Wave  Energy  Growth  in  Shallow  Water 

In  Chen  and  Wang's  (1983)  model,  the  generation  term  is  again  defined  by 
Eq.(2.47).  In  order  to  be  consistent  with  the  deepwater  model,  we  propose  to  use 
the  same  functional  form  here  for  local  wave  generations  as  for  deepwater  case: 


G,n  -  { 


0.1  s f{^y[E,{f)- cos' {d  -  e^)  -  E{f,9)  ],     if  1^  -  e^\<  ir/2 

^  ^  (3.22) 

0,  elsewhere 

where  E,{f)  is  the  duration  unlimited  shallow  water  spectrum.   The  spectral  evo- 
lution based  upon  Equation  (3.22)  at  a  duration-limited  Ccise  can  be  expressed  as 

E'^'ifJ)   =  El^,{f,e)exv{-0.1sf{^rAt] 

+  £,{/){!-  exp[-0.1sf{^)*At]}-cos\9  -  e^)         (3.23) 
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This  equation  is  actually  used  in  the  UCWP  shallow  water  model,  and  it  is  used  to 
compute  the  wave  components  for  UiqJC  >  0.9  such  that  it  will  not  contajninate 
swells. 

Currently,  our  knowledge  on  E,{f)  is  limited.  In  the  present  model  we  adopted 
the  TMA  spectral  form  proposed  by  Bouws  et  al.(1983): 

na^  ^  exp ( r-) 

^'^^^  =  (2^"^f-4(^-/'^^'l  ^       ^     2    a^  /^  ^  ^$(fc,;,)  (3.24) 

where  $(A;,/i)  is  a  transfer  fimction  defined  earlier  in  Eq.(2.49): 

_  fc^(/,oo)C,(/,oo) 

The  spectral  parameter  a  is  estimated  in  the  model  by  the  following  formula: 

a  =  0.0156  (EkhlLf''^  (3.25) 

where  A:^  is  the  corresponding  wave  number  at  peak  spectral  frequency,  /^,  which  in 
turn  is  related  to  the  resonant  frequency  as  fm  =  0.88  /q.  Equation  (3.25)  describes 
a  best-fit  curve  of  the  upper  boim.d  values  of  a  computed,  beised  on  field  data,  by 
Vincent  (1984).  The  spectral  parameters  7,  and  aat  are  also  adopted  from  Vincent 
(1984): 

^  =  2.47(^^^)°^  (3.26) 

'  0.07     if  /  <  /^, 
Oai  =  {  (3.27) 

,  0.09     if  /  >  /^. 

However,  since  a  significant  amount  of  wave  energy  is  continuously  being  dissipated 
in  the  low  frequency  range  under  the  influence  of  shallow  water,  the  energy  densities 
in  this  rEOige  may  not  be  sufficiently  lasge  to  maintain  the  spectral  peak  near  the 
resonant  frequency.  Vincent  (1984)  suggested  that  the  peak  spectral  frequency 
should  have  a  lower  bound 


....  >  ^('-^" 

where  A  lies  between  0.8  and  1.0 


/m  >   -(J)'/'  (3.28) 
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3.3.2     Wave  Energy  Dissipation  in  Shallow  Water 

The  dissipation  of  wave  energy  in  the  shallow  water  may  be  attributed  to  the 
following  effects:  (i)  flow  turbulence,  (ii)  bottom  friction,  (iii)  bottom  percolation, 
(iv)  breaking  due  to  depth  limitation  and  (v)  breaking  due  to  energy  saturation  of 
spectral  components. 

The  dissipation  mechainism  due  to  flow  turbulence  is  essentially  the  same  as 
in  deep  water,  or  Eqs.(3.17)  and  (3.18)  still  apply.  The  bottom  friction  eff'ect  is 
estimated  from  Eqs.(2.50)  and  (2.51).  The  dissipation  of  wave  energy  due  to  perco- 
lation of  porous  sand  bottom  is  estimated  from  Eq.(2.53).  However,  this  quantity 
is  found  to  be  very  small  and  can  be  neglected.  The  water  depth  limited  condition 
is  given  by  Eq.(2.57). 

The  limiting  of  spectral  growth  due  to  energy  saturation  is  based  on  Eq.(2.56), 
or 

A{f)  =  Bg'{2^f)-'^{k,h] 

where  B  is  chosen  equal  to  0.073  in  this  study. 

3.4     Wind  Model 

Wind  is  aji  important  factor  for  generating  water  siirface  waves.  The  wind 
information  required  as  input  to  a  wave  model  is  the  time  series  of  the  overwater 
wind  vector  fields.  The  wind  information  can  be  obtained  either  from  the  field 
mccisurement  or  from  a  wind  model.  Since  measurement  of  the  winds  over  aji  entire 
water  surface  region  is  almost  impossible,  numerical  generation  of  the  wind  field  is 
generally  the  only  means  to  prepare  the  wind  information  as  input  to  a  wave  model. 

The  numerical  model  utilized  to  generate  the  wind  information  here  follows  the 
method  described  by  Clarke  cind  Hess  (1974).  It  generates  wind  field  based  upon 
the  input  of  sea  surface  pressure  distribution.  The  proceduces  are  summarized  as 
follows: 
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(i)  The  geostrophic  gradient  winds  are  produced  directly  from  the  pressure  gradient 
information  by  the  following  equation  (Byers,  1974) 

uf+  — = -^  3.29 

T  Pa  dn 

where  u=|u|  is  the  magnitude  of  wind  speed,  r  and  n  are,  respectively,  the  radius  of 
curvature  and  the  directional  normal  of  the  isobars,  Pa  is  air  density,  f=  211  sin  (^  is 
the  local  component  of  the  planetary  vorticity  normal  to  the  earth's  surface,  or  the 
Coriolis  parameter  with  Cl  the  angular  velocity  of  earth  rotation  and  <j>  the  latitude, 
cind  p  is  the  atmospheric  pressure.  If  the  geostrophic  gradient  wind  is  to  be  solved 
in  the  Cajtesian  coordinates,  (x,  y),  dp/dn  can  be  obtained  from  the  expression 

an  ox  ay 

and  the  radius  of  curvature,  r,  is  obtained  from 

,,d(p.         .  ,     d(p       dtp  dx       d<p  dy       dp         „        dp 
r  =  l/K7  >    with  -^  =  -^— +  -^-^  =  -^cos5, +  -!^sin^, 
di  dt        dx  di       dy  di        dx  dy  " 

where  £  is  the  distance  along  the  path  where  air  flows,  and  Og  is  the  downwind 

direction  tangent  to  the  isobar.  It  should  be  noted  that  in  the  Northern  Hemisphere 

a  positive  r  indicates  the  cyclones  and  a  negative  r  indicates  the  anticyclones.  The 

geostrophic  gradient  winds  solved  from  Eq.(3.29)  are  then  used  to  compute  the 

winds  at  a  reference  level  above  the  sea  surface. 

(ii)     The  surface  winds  are  then  obtained  by  the  corrections  of  frictional  effects. 

The  wind  speed,  U^,  and  wind  direction,  6^,  are  adjusted  according  to  the  following 

equations: 

^  =  —  In  -  (3.30) 

U,       0.4       20  ' 


and 


d^  -  6g  =  s\n-\Bol —^)  (3.31) 
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where  the  frictional  velocity,  U.,  and  the  frictional  height,  Zq,  are  computed  by 
(Clarke  and  Hess,  1974) 

i^^y  -  {BoV  +  M^)  -  Ao]'  (3.32) 

U,  I  Zq 

and  Eq.(2.42),  or 

0.1525  Uj 

Zq  =  —- h0.0144—^ -0.00371,       (in  cgs  units) 

U,  g 

where  Aq  ajid  Bq  are  non-dimensional  parameters,  depending  on  horizontal  temper- 
ature gradient  and  vertical  heat  flux  in  the  surface  boundary  layer.  At  adiabatic 
condition  Aq  and  Bq  are  contants  and  equal  to  0.8  and  3.5,  respectively.  More  in- 
formation about  the  values  Aq  and  Bq  in  the  nonadiabatic  condition  can  be  found 
in  Resio  et  al.(l982),  and  Large  and  Pond  (1982). 


CHAPTER  4 
MODEL  PERFORMANCE 


4.1     Introduction 


The  new  UCWP  model  is  to  be  tested  here  for  hypothetical  and  real  cases. 
The  hypothetical  cases  consider  the  wave  generation  at  either  fetch  limitation  or 
duration  limitation  imder  stationary  wind  fields.  Only  the  deepwater  submodel 
is  tested  for  the  hypothetical  cases.  The  real  cases  consider  the  hindcastings  of 
deepwater  waves  in  the  Northwest  Atlaintic  Ocean  and  shallow  water  waves  along 
the  East  Florida  Coast.  Both  the  deepwater  and  shallow  water  submodels  axe  tested 
in  the  real  cases. 

In  Section  4.2  the  UCWP  deepwater  submodel  was  tested  against  those  of  the 
WIS  model,  the  SMB  model,  and  Hasselmann's  (1976)  parametric  model.  Such  tests 
were  made  for  either  fetch-limited  or  duration-limited  condition  at  selected  constajit 
wind  speeds.  In  Section  4.3  the  UCWP  model  Weis  tested  against  the  measured  data 
from  NOAA  buoy  deployed  in  the  Northwest  Altantic  Ocean  and  from  CDN  gages 
installed  along  the  East  Florida  Coast.  In  Section  4.4  the  relationship  between  wind 
and  wave  for  a  fully  aroiised  sea  is  studied.  In  Section  4.5  the  relative  importance  of 
each  factor  affecting  deep  and  shallow  water  spectral  transformations  is  analyzed. 
Finally,  in  Section  4.6  various  values  of  bottom  friction  coefficient  were  tested  in 
the  UCWP  shallow  water  model  to  examine  their  effect  on  the  wave  transformation 
process. 
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4.2     Model  Performance  in  Hypothetical  Cases 

In  this  section,  the  UCWP  model  is  compared  with  three  standard  wave  hind- 
casting  models,  the  WIS  model,  the  SMB  model,  and  the  Hasselmann's  (1976) 
pajajnetric  model,  for  fetch-limited  ajid  duration-limited  deepwater  wave  genera- 
tions. The  deepwater  WIS  model  and  Hasselmann's  parametric  model  have  been 
reviewed  in  Chapter  2.  The  SMB  model  is  known  as  a  semiempirical  method  for 
deepwater  wave  prediction;  it  was  first  developed  by  Sverdrup  and  Munk  at  1947 
cind  then  revised  by  Bretschneider  at  1952  and  1958  with  additional  empirical  data. 
The  wave  properties  used  for  comparisons  are  the  significant  wave  height  and  modal 
period. 

The  UCWP  deepwater  submodel  wcis  tested  first  in  fetch-limited  cases  with 
constant  wind  speeds  equal  to  20,  30,  40,  and  50  knots,  respectively.  The  results  of 
significant  wave  height  H,  and  modal  period  Tm  are  plotted  against  fetch  at  20,  50, 
100,  150,  200,  250,  300,  350,  400,  450,  500,  and  550  kilometers  in  Fig.  4.1. 

According  to  the  SMB  model  the  significant  wave  height,  H,,  and  modal  period, 
r^  along  a  fetch  F  axe  estimated  from  the  following  formulas,  respectively: 

^=  0.283  tanh[0.0125(^)^/n  (4.1) 


y^=7.54tanh[0.077(^)^/'']  (4.2) 

Cio  t/io 

where  11^,  the  adjusted  wind  speed,  is  defined  as  a  power  function  of  Uxq  (Shore 
Protection  Manual,  1984): 

Ua  =  0.71  C/io  (iri  meter-second  units)  (4.3) 

According  to  Hasselmarm's  parametric  model  the  significant  waves  and  modal 
periods,  along  a  fetch  F,  are  estimated  from  the  following  equations: 

^=0.0016(|^)"'  (4.4) 
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^=  0.286  (^)^/^  (4.5) 

In  the  classical  SMB  and  parametric  models,  the  adjusted  wind  speed  Ua  is 
equal  to  the  wind  speed  Uiq.  However,  it  is  found  later  that  a  correction  to  Ua,  is 
necessary  since  both  the  classical  SMB  model  and  the  Hasselmann's  model  tend  to 
imderestimate  the  wave  heights.  The  utilization  of  Ua  according  to  Eq.(4.3)  in  the 
both  models  is  suggested  by  the  Corps  of  Engineers,  U.S.  Army,  and  Ua  obtained 
from  the  equation  will  be  greater  than  Uiq  if  Uw  is  greater  than  4  meters  per  second. 
Therefore,  by  using  the  adjusted  wind  speed  Ua  as  defined  in  Eq.(4.3),  the  SMB 
and  parametric  models  will  predict  higher  wave  heights  than  the  classical  ones  if 
Uio  is  greater  than  4  meters  per  second. 

The  significant  wave  heights  and  modal  periods  obtained  from  the  WIS  model, 
the  SMB  model,  and  Hasselmann's  parametric  model  are  compared  to  the  UCWP 
model  in  Figure  4.1.  At  low  and  moderate  wind  speeds  (20,  30  and  40  knots) 
and  large  fetch  (greater  than  400  km),  the  significant  wave  heights  obtained  from 
the  UCWP  model  are  quite  consistent  with  these  obtained  from  the  WIS  model. 
At  high  wind  speeds  (40  and  50  knots)  and  small  fetch  (less  than  400  km),  the 
significant  wave  heights  obtained  from  the  UCWP  model  are  the  highest  among  the 
four  models,  while  the  significant  heights  obtained  from  the  WIS  model  and  the 
parametric  model  agree  with  each  others.  As  fetch  increases,  the  significant  wave 
heights  obtained  from  the  UCWP  model,  WIS  model,  and  SMB  model  tend  to  reach 
to  a  height  limit  due  to  energy  saturation.  The  parametric  model,  however,  is  not 
bounded. 

The  modal  periods  obtained  from  the  UCWP  model  au-e  consistent  with  those 
obtained  from  the  SMB  model  at  low  and  moderate  wind  speeds  (20  ajid  30  knots). 
At  high  wind  speeds  (40  and  50  knots),  the  modal  periods  obtained  from  the  UCWP 
model  are  much  smaller  than  the  other  three  models.  For  high  wind  speeds  (40 
and  50  knots)  zrnd  small  fetch  (less  than  400  km),  the  modal  periods  computed 
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from  the  WIS  model  and  the  parametric  model  agree  with  each  other.  At  large 
fetch,  the  modal  periods  obtained  from  the  UCWP  model,  the  WIS  model,  and  the 
SMB  model  again  tend  to  reach  an  upper  limit  due  to  saturated  conditions.  The 
pzwametric  model  does  not  have  aji  upper  bound  on  modal  period. 

The  UCWP  deepwater  submodel  was  compared  next  in  the  duration-limited 
cases  with  the  WIS  model,  the  SMB  model,  and  Hasselmann's  model  at  the  constant 
wind  speed  of  20,  30,  40,  and  50  knots,  respectively.  The  comparisons  are  shown  in 
Figiire  4.2. 

The  significant  wave  height  and  modal  period,  respectively,  for  a  duration  t  axe 
estimated  in  the  SMB  model  from  the  following  formulas: 

^—^  =  0.283  tanh[0.00087(^)^/«]  (4.6) 

^  =  7.54  tanh[  0.016(^)^/«  ]  (4.7) 

where  the  adjusted  wind  speed,  Ua,  is  defined  in  Eq.(4.3). 

The  significant  waves  and  modal  periods  for  a  duration  t,  are  estimated  in  the 
Hasselmann's  parametric  model  from  the  following  equations: 

^  =  0.000082  (^)^/^  (4.8) 

^=0.06(^)^/^  (4.9) 

It  is  seen  in  Figure  4.2  that  at  short  duration  (less  than  20  hrs)  consistent  sig- 
nificant wave  heights  are  found  between  the  UCWP  and  SMB  models,  and  between 
the  WIS  and  Hasselmann's  models.  At  high  wind  speeds  (40  and  50  knots)  and 
short  duration  (less  than  20  hrs),  the  significant  wave  heights  obtained  from  the 
UCWP  and  SMB  models  are  found  to  be  much  greater  than  the  WIS  and  para- 
metric models.  As  a  matter  of  fact,  the  significant  wave  heights  obtained  from  the 
UCWP  model  are  the  highest  among  the  four  models  at  short  duration.  At  long 
duration  (greater  than  20  hrs),  both  the  SMB  and  parametric  models  predict  much 
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higher  waves  than  the  UCWP  and  WIS  models.  At  long  duration,  wave  heights 
generated  from  the  UCWP  and  SMB  models  have  upper  limits.  The  significant 
heights  generated  from  the  UCWP  model  show  'overshoot'  patterns  in  Fig.  4.2.  In 
other  words,  waves  become  over  saturated  before  settle  to  a  stable  saturation  con- 
dition. The  largest  'overshoot'  0.2  m  associated  at  30  knots  winds;  it  is  still  very 
small  (only  3  %)  compared  with  a  6.6  m  height  at  saturation.  The  'overshoot'  of 
wave  height  is  overall  small  but  it  is  essential  to  gain  enough  energy  for  establishing 
the  saturation  condition  in  the  UCWP  model.  The  'overshoot'  does  not  occiir  in 
the  WIS  model  nor  in  the  SMB  and  parametric  models.  None  of  the  four  models 
shows  overshoot  behavior  in  the  fetch-limited  cases. 

The  significant  wave  heights  predicted  from  the  WIS  model  at  40  knots  wind 
show  a  peculiar  periodic,  irregular  damping  pattern  when  duration  becomes  longer 
thaji  40  hrs.  This  damped  wave  height  which  is  unlikely  to  occurr  in  the  real 
situation  ranges  from  10.6  m  to  13.4  m  and  repeats  its  pattern  every  12  hrs.  It 
appears  that  it  is  caused  by  an  improper  choice  of  saturation  spectrum  E,[f),  or 
Eq.(2.32),  which  function  is  to  limit  the  growth  of  waves  in  the  WIS  model. 

The  modal  periods  obtained  from  the  four  models  do  not  differ  appreciably  to 
each  others  for  small  duration  (less  than  10  hrs).  For  moderate  duration  (between 
10  hrs  and  40  hrs),  the  modal  periods  obtained  from  the  WIS  model  agree  with  the 
SMB  model  at  low  wind  speed  (20  and  30  knots)  but  become  closer  to  that  predicted 
by  the  parametric  model  at  high  wind  speed  (40  and  50  knots).  For  long  duration 
(greater  than  40  hrs),  consistent  modal  periods  are  found  between  the  WIS  and 
SMB  models.  At  long  duration,  the  modal  periods  obtained  from  the  parajnetric 
model  are  the  largest  among  the  four  models.  At  long  duration  and  high  winds  (40 
aind  50  knots),  the  modal  periods  obtained  from  the  UCWP  model  are  found  to  be 
much  less  than  the  other  three  models. 
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Figure  4.3:  Shallow  water  grid  system. 


4.3     Model  Performance  in  Real  Cages 


The  UCWP  model  was  first  tested  against  the  bench  mark  case  of  the  period 
December  21,  1983,  to  January,  1984  to  offer  comparisons  with  the  WIS  model 
and  field  data.  For  the  deepwater  model,  the  test  conditions  including  grid  size, 
time  step,  etc.,  are  identical  to  that  described  in  the  WIS  model  testing  as  shown 
in  Section  2.7.  For  the  shallow  water  model  the  grid  system  is  given  in  Fig.  4.3. 
The  grid  system  covers  an  area  300  kilometers  long  -  from  mid  Georgia  to  the 
northern  edge  of  Bahsima  Islands,  and  about  100  kilometers  wide  -  from  shoreline 
to  approximately  600  m  contour  line  offshore. 

The  wind  information  was  generated  by  the  wind  model  with  daily  synoptic 
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surface  pressure  as  input.  Since  the  surface  chart  is  given  daily  at  7:00  a.m.,  E.S.T. 
only,  linear  extrapolations  are  used  for  winds  at  other  times.  The  computed  wind 
directions  and  speeds  for  January  1984  at  grid  point  30.0°N  and  77.5''W  are  shown 
in  Fig.  4.4  together  with  the  measured  winds  from  NOAA  buoy  #41006,  located 
at  29.3''N  and  77.3°W.  The  correlation,  6,  given  in  Fig.  4.4  has  been  defined  in 
Eq.(2.58).  The  root  mean  square  difference  in  wind  direction,  Ormse-,  given  in 
Fig.  4.4  is  defined  as  follows: 

Ormse  =  {T,li^..c)i  -  [e.M'V^VN  (4.10) 

where  d^^c  is  the  computed  wind  direction,  9,ju^m  is  the  measured  wind  direction,  A'^ 
is  the  sample  size,  and  the  overhead  bar  indicates  the  operation  of  arithmatic  mean. 

The  wind  strength  exhibits  a  periodical  behavior  with  a  periodicity  of  roughly 
equal  to  7  days.  The  wind  direction  associated  with  high  winds  are  found  to  be 
more  stable  and  to  have  better  agreement  with  the  predicted  values.  The  small 
winds  are  less  stable,  thus,  less  predictable. 

The  measured  wind  speeds  are  seen  in  Fig.  4.4  to  be  slightly  scattered  around 
the  computed  ones  within  ±10  knots.  This  is  because  the  measured  winds  are 
reported  in  an  every  hour  schedule  while  the  computed  winds  are  from  the  wind 
model  at  an  every  24  hrs  interval.  The  fluctuations  of  either  magnitude  or  direction 
of  hourly  recorded  winds,  as  seen  in  Fig.  4.4,  may  be  caused  by  abrupt  changes  of 
temperature  between  the  air  and  sea  water,  and  between  the  stratified  air  layers 
above  the  sea  surface.  For  simplicity  the  temperature  effect  is  not  include  in  the 
wind  model  which  assumes  a  neutral  condition. 

The  wind  information  used  as  input  to  the  wave  model  are  the  computed  winds 
adjusted  by  the  meaisured  values  at  the  buoy.  That  is,  the  wind  speed  at  any  grid 
point  is  the  computed  value  multiplied  by  a  factor  which  is  the  ratio  of  measured 
and  computed  ones.  However,  this  factor  is  limited  in  the  range  between  0.5  and 
1.5  to  limit  the  influence  of  abnormal  data  values. 
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Figure  4.4:    Comparison  of  hindcasted  and  measured  winds  at  the  Buoy  #41006 
location. 
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Figure  4.5  shows  the  comparison  of  the  time  series  of  the  wave  energy  spectral 
evolution  process  of  the  entire  month.  Figure  4.6  shows  the  comparisons  between 
the  meaisured  and  computed  modal  periods  and  significant  wave  heights  at  the  buoy 
location.  Figure  4.7  shows  the  comparison  between  the  measured  and  computed 
energy  scale  parameters  at  the  buoy  location.  The  model  performs  well  for  this 
data  set.  The  good  agreement,  however,  should  be  viewed  as  partially  ajtificial  as 
the  model  was  calibrated  on  this  data  set.  The  model  outputs  in  the  ensuing  cases 
axe  data  independent. 

In  the  shallow  water  model,  only  one  empirical  value,  the  bottom  friction  coef- 
ficient, needs  to  be  defined  as  input.  Shemdin  et  al.(1978)  suggested  values  ranging 
from  0.006  to  0.05  along  the  east  Florida  coast.  A  medium  value  0.02  was  chosen 
for  the  model.  Directional  wave  spectra  from  deepwater  model  are  used  as  input 
boundary  conditions.  The  wind  field  is  spatially  linearly  extrapolated  from  the  ad- 
justed output  from  the  wind  model  as  discussed  earlier.  All  the  input  conditions  are 
upgraded  at  2  hrs  interval  but  the  time  step  in  the  shallow  water  model  is  chosen 
as  24  minutes  to  insure  numerical  stability.  The  dissipation  of  wave  energy  due  to 
percolation  is  not  considered  in  the  shallow  water  model  since  it  was  found  to  be  a 
minor  effect  to  the  wave  transformation  process  in  the  water  region  neax  the  Eeist 
Florida  Coast  (Shemdin  et  al.,  1978). 

For  the  case  of  January,  1984,  the  hindcast  results  from  the  shallow  water  model 
eire  compared  with  the  directional  wave  information  collected  at  the  offshore  Kings 
Bay  Station  (St. #5  in  Fig.  2.4).  Figure  4.8  shows  the  comparison  of  the  time  series 
of  the  one-dimensional  spectra.  By  comparing  Fig.  4.8  with  Fig.  4.5,  the  measured 
deepwater  and  shallow  water  wave  spectra,  we  noticed  the  significant  reduction  in 
wave  energy  and  the  reduction  is  particularly  severe  in  low  frequency  components. 
This  point  will  be  further  discussed  later.  The  time  series  of  the  marginal  directional 
spectra,  which  are  obtained  by  integrating  the  two-dimensional  spectra  over  the 
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Figure  4.5:    Comparison  of  NOAA  and  UCWP  wave  spectra  for  January  1984  at 
the  Buoy  #41006  location. 
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Figure  4.6:  Comparisons  of  NOAA  and  UCWP  modal  periods  and  significant  wave 
heights  at  the  Buoy  #41006  location. 


82 


10- 


flH, 


(UCWP) 


u.  i<: 

+    2n>H3 

6 

= 

0.44 

A    3H>Hsi2H 

6 

= 

0.43 

♦           Hs 

i3M 

5 

= 

0.67 

0.  10 

A 

«          A 
O        A 

+ 

+ 

0.08 

+  +                ^ 

+ 

A 

^                 + 

^  ^              ^  «^ 

+ 

^ 

♦ 

A 

*  A                  -M-                              ^ 

0.06 

A      A    aO 

A 

A 

^     '^   ^ 
>    A            <> 

A 

0.04 

A*  *-^.  ^ 

6a   * 

0.02 

A        -^    *+    0      -^       * 

^    aI   +^ 

+ 

A+ 

4:a  + 

^#+*A^                      AAA       V 

4-                       A 
^A 

A 

0.0 

1                                      1 

1 

1 

0.0 


0.02 


10 


0.04 


0.06 


0.08 


0.  1 


,  (NOAA) 


Figure  4.7:  Comparison  of  NOAA  and  UCWP  energy  scale  parameters  at  the  Buoy 
#41006  location. 
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Figure  4.8:  Comparison  of  CDN  and  UCWP  wave  spectra  for  January  1984  at  the 
Kings  Bay  gage  location. 
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frequency  domain,  cire  shown  and  compared  in  Fig.  4.9.  Finally,  the  measured  and 
computed  mean  wave  directions,  modal  periods,  and  significant  wave  heights  are 
summarized  in  Fig.  4.10. 

In  the  case  illustrated  above,  the  winds  were  predominantly  induced  by  high 
pressure  systems  which  are  large  scale  and  relatively  stationary.  In  the  following 
case,  wave  hindcasting  was  performed  in  September  auid  October,  1984  to  assess 
the  model  performance  for  low  pressure  system  and  other  coastal  locations. 

In  September  and  October,  1984,  three  hurricanes  occurred  in  the  Mid  Atlantic 
Ocean  (see  Fig.  4.11).  The  three  hurricanes  are  Hurricane  Diana,  Isidore,  and 
Josephine.  The  cyclones  found  during  these  hurricanes  are  good  examples  of  low 
pressure  system  found  in  September  ajid  October;  they  are  small  scale  and  relatively 
nonstationajy.  Synoptic  charts  showing  low  pressure  systems  of  Hurricane  Diana, 
Isidore,  and  Josephine  are  given  in  Figs.  4.12(a),  4.12(b),  and  4.12(c),  respectively. 
All  three  cylcones  generate  strong  northeast  winds  and,  hence,  large  waves  heading 
to  the  southeastern  coast  of  the  United  States. 

Hurricane  Diana  came  near  the  Florida  east  coast  in  the  morning  of  Septem- 
ber 8,  1984.  After  looped  off  the  northeast  Florida  coast  for  almost  36  hours,  it 
moved  north  toward  the  coast  of  North  Carolina.  The  highest  wave  recorded  at 
the  Marineland  CDN  station  during  the  visit  of  Hurricane  Diana  reached  3  m  in 
terms  of  the  significant  wave  height.  Hurricane  Isidore  came  ashore  near  Miami 
Beach,  Florida,  during  the  morning  of  September  27,  1984.  It  then  moved  north- 
ward ax;ross  central  Florida  and  returned  to  the  Atlantic  Ocean  near  the  border 
between  Florida  and  Georgia.  The  highest  wave  caused  by  the  Hurricane  winds 
was  recorded  at  the  CDN  wave  gauge  station  at  Cape  Canaveral  and  reached  3.2 
m  in  terms  of  the  significant  waves.  Hurricane  Josephine  v/as  developed  in  the  Mid 
Atlajitic  Ocean  on  October  10,  1984.  The  path  of  the  hurricane  is  nearly  parallel 
to  the  east  coast  of  North  America.   The  closest  distance  between  the  hurricane's 
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Figure  4.9:    Comparison  of  CDN  and  UCWP  marginal  directional  spectra  at  the 
Kings  Bay  gage  location. 
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Figure  4.10:  Comparisons  of  CDN  and  UCWP  average  wave  direction,  modal  peri- 
ods and  significant  wave  heights  at  the  Kings  Bay  station. 
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Figure  4.11:  Tracks  of  Hurricanes  Diana,  Isidore,  and  Josephine  (1984). 
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Figure  4.12:  (a)  Weather  map  showing  northeast  winds  generated  by  the  Hurricane 
Diana  (9/8/84). 
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Figure  4.12:  (b)  Weather  map  showing  northeast  winds  generated  by  the  Hurricane 
Isidore  (9/27/84). 
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Figure  4.12:  (c)  Weather  map  showing  northeast  winds  generated  by  the  Hurricane 
Josephine  (10/11/84). 


91 

center  and  the  east  coast  of  Florida  is  about  800  kilometers.  The  hurricane  never 
threatened  Florida  coast.  However,  due  to  its  track,  strong  northeasters  were  expe- 
rienced along  the  coast  of  Florida  and  large  waves  were  observed.  The  significant 
wave  height  recorded  2.2  m  at  the  Maxinelcind  CDN  station. 

Two  large  scale  high  pressure  weather  systems  were  also  developed  over  the 
Northwest  Atlantic  Oceaji  in  September  and  October  1984.  The  first  one  was  orig- 
inally formed  in  the  central  area  of  North  America  ajid  later  moved  eastward  into 
the  Atlantic  Ocean  on  the  22nd  of  September.  It  stayed  over  the  water  regions  near 
the  coast  of  North  Cajolina  for  about  three  days  while  generating  persistent  strong 
winds  blowing  to  the  waters  surrounding  the  southern  coast  of  North  America.  Dur- 
ing the  presence  of  the  ainticyclones  the  highest  significajit  wave  height  wcls  recorded 
at  the  Kings  Bay  Station  at  2.5  meters.  The  second  large  scale  high  pressure  system 
appeared  over  the  coastal  waters  neax  the  coast  of  North  Carolina  on  the  25th  of 
October.  It  lingered  there  for  about  two  days  before  moving  northward.  During 
the  anticyclone  period  the  significant  wave  height  reached  2.1  m  at  the  Marinelajid 
CDN  station. 

The  comparisons  of  hindcasted  values  ajid  measured  data  are  given  in  a  series 
of  graphics.  Figure  4.13  shows  the  comparisons  between  modal  periods  aind  signif- 
icajit wave  heights  at  the  deepwater  buoy  location.  Figiires  4.14(a),  4.14(b),  and 
4.14(c)  shows  the  comparisons  between  average  wave  directions,  modal  periods, 
ajid  significajit  wave  heights  at  the  shallow  water  stations  of  Kings  Bay,  Marineland 
ajid  Cape  Canaveral,  respectively.  A  summary  of  compcirisons  of  the  computed 
and  measured  unidirectional  frequency  spectra  at  the  buoy  #41006  location,  for  all 
Januajy,  September  asid  October  1984  data  is  given  in  Appendix  C. 

Table  4.1  summarizes  the  correlations  between  the  computed  and  mccisured  val- 
ues for  significant  wave  height  and  modal  period  at  the  buoy  #41006  location.  From 
the  table  it  is  evident  that  higher  correlations  are  achieved  for  both  wave  height 
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Figure  4.13:  Comparisons  of  NOAA  and  UCWP  modal  periods  and  significant  wave 
heights  at  the  Buoy  #41006  location. 


Table  4.1:  Summary  of  Correlations  between  computed  and  meaisured  values  of  H, 
and  Tm  at  the  buoy  #41006  location. 


station 

month 

WIS  model 

NOAA 

buoy 
#41006 

Jan. (84) 

"Hm,c  Ht,m 

Tm,c  Tm,fa 

0.49 

0.49* 

-0.22 

-0.21* 

Sep.-  Oct. (84) 

0.88 

0.90* 

0.30 

0.50' 

station 

month 

UCWP  model 

NOAA 

buoy 

#41006 

Jan. (84) 

Ht.c  Ht,m 

Tm,c  Tm,,m 

0.85 

0.86* 

0.38 

0.53- 

Sep.-  Oct. (84) 

0.90 

0.92' 

0.14 

0.51' 

*swell  data  {Hi,c<2m,Tm.,m>8sec)  are  excluded  in  correlation  calculations 
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Figure  4.14:   (a)  Comparisons  of  CDN  and  UCWP  average  wave  direction,  modal 
periods  and  significant  wave  heights  at  the  Kings  Bay  gage  location. 
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Figure  4.14:   (b)  Comparisons  of  CDN  and  UCWP  average  wave  direction,  modal 
periods  and  significant  wave  heights  at  the  Maxineland  gage  location. 
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Figure  4.14:   (c)  Comparisons  of  CDN  and  UCWP  average  wave  direction,  modal 
periods  and  significant  wave  heights  at  the  Cape  Canaveral  Station. 
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Table  4.2:  Summary  of  Correlations  between  computed  and  measured  values  of  H, 
and  r^  at  three  shallow  water  stations. 


station 

month 

UCWP 

model 

Kings  Bay 

Jan.(84) 

^H..c 

H  9,m. 

^T,n.c 

Tvn.,m, 

0.77 

o.ir 

0.45 

0.34' 

Sep. 

Oct. (84) 

0.81 

0.79* 

-0.05 

0.13* 

Marineland 

0.92 

0.93* 

0.16 

0.57* 

Cape  Canaveral 

0.77 

0.80' 

0.37 

0.22* 

*swell  data  [Hs,c<im.,Tm.,m.>\Qsec)  are  excluded  in  correlation  calculations 


and  wave  period  when  swell  data  are  excluded.  The  comparisons  of  correlations  also 
show  that  the  UCWP  deepwater  submodel  yields  higher  and  more  consistent  cor- 
relations than  the  WIS  deepwater  model.  Table  4.2  shows  the  correlations  between 
the  computed  and  measured  signficant  wave  heights,  as  well  as  modal  periods,  at  the 
shallow  water  stations  of  Kings  Bay,  Marineland  and  Cape  Canaveral.  The  correla- 
tions between  the  computed  and  measured  significant  wave  heights  do  not  change 
appreciably  whether  or  not  swell  data  are  included  in  the  computation  of  corre- 
lations. However,  the  correlations  of  the  computed  and  mecisured  modal  periods 
improved  significantly  if  swell  data  are  excluded  in  the  correlation  calculations. 

In  summary,  the  model  performed  well  in  wave  height  and  direction  hindcastings 
especially  for  severe  weather  cases.  Under  no  wind  or  weak  wind  condition,  the 
model  tended  to  underestimate  the  waves.  The  overall  performance  on  modal  period 
hindcasting  was,  at  best,  fair.  The  computed  modal  periods  are  consistent  with  the 
measured  ones  for  large  waves  but  fare  poorly  for  small  waves.  Also,  for  mixed 
swell  and  sea  conditions,  modal  period,  Tm,  appears  to  be  a  poor  indicator  of  wave 
property. 

4.4     Wind- Wave  Relationship  in  a  Fully  Developed  Sea 

In  Section  4.2  the  significant  wave  heights  and  modal  periods  computed  from 
the  UCWP  model  are  shown  to  converge  to  their  respective  saturation  limits  at 


97 

different  wind  speeds  in  both  fetch  and  duration  dependent  cases.  Based  upon 
these  wave  data,  Fig.  4.15  shows  on  log-log  scale  the  plots  of  the  significant  wave 
height  versus  wind  speed,  and  the  modal  period  versus  wind  speed  for  different 
fetchs  and  durations.  It  is  seen  in  the  plot  of  H,  versus  Uiq  both  fetch-limited  and 
duration-limited  curves  tend  to  approach  to  a  single  limit  for  large  fetch  and  long 
duration,  which  indicates  a  fully  developed  condition  for  wind  generated  waves.  For 
duration  dependent  case,  it  is  interesting  to  observe  here  that  wave  height  and  wave 
period  approach  to  their  respective  upper  limits  at  approximately  the  same  rate.  If 
the  fetch  is  unlimited,  the  wave  can  be  considered  virtually  saturated  (both  height 
and  length  reached  limit)  when  duration  exceeds  60  hours.  For  fetch  dependent 
case,  on  the  other  hand,  wave  height  reaches  its  upper  limit  before  the  wave  period, 
that  is,  as  the  fetch  increases  beyond  600  km,  wave  height  virtually  reaches  its  limit 
and  becomes  more  difficult  to  grow  while  wave  period  can  still  increase  easily  to 
some  extent  utill  it  reaches  the  upper  limit. 

In  fetch  dependent  case,  an  empirical  equation  for  H,  similar  to  that  of  SMB 
formula  is  proposed  here: 

7S^=<^i[7^-^i-tanh(c3|^)]-  (4.11) 

where  ci,  cz,  c^  are  constants,  F^  is  the  fetch  corresponding  to  maximum  overshoot 
(see  Fig.  2.4),  and  Ri  is  a  overshoot  enhancement  factor.  The  overshoot  enhance- 
ment parameter  7i  and  the  width  parameter  Oi  are  functions  of  Uiq.  The  above 
equation  has  similar  functional  form  8is  the  parametric  model  when  F I F,  is  small, 
i.e., 

where  xi  's  a  proportional  constant.  For  large  F/Fg,  Ri  approaches  to  unity  and 
Eq.(4.11)  can  be  expressed  as 

7^  =  ^1(7^)".     forF»F,  (4.12) 
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Figure  4.15:  Plots  of  (a)  H,  versus  Uw,  and  (b)  T^  versus  Uw  at  various  fetch  and 
duration. 
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Thus,  Fa  can  be  treated  as  the  practical  upper  limit  of  the  fetch  influence. 

In  duration  dependent  Ccise,  an  empirical  equation  for  H3  may  be  expressed  as 

f^  =  C4[^-i22-tanh(c5f)p  (4.13) 

where  C4,  C5,  cq  are  constants,  t,  is  the  duration  required  for  maximum  overshoot, 
and  R2  is  again  a  overshoot  enhancement  factor.  Similaxly,  when  t  <^  t,  Eq.(4.13) 
can  be  approximated  by  the  expression: 

^-^  =  X2[j^Y\     fori<«, 

where  X2  is  a  proportional  constant,  and  when  f^t,  Eq.(4.11)  can  be  expressed  as 

7^  =  '^4(^r,     forf»f,  (4.14) 

At  the  fully  developed  condition,  for  fetch  and  duration  unlimited  case  Eq.(4.12) 
shall  be  identical  to  Eq.(4.14),  or 

93a,<K 


u}^ 


criU^oY'  (4.15) 
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where  H,oo  indicates  the  significant  height  at  the  fully  developed  sea.  If  this  fully 
developed  condition  is  approximated  by  a  upper  bound  of  H,  —  Uio  curves,  drawn  as 
a  solid  line  in  Fig.  4.15,  we  find  cj  «  0.738  (m/sec)^/^  and  cg  «  -1/3,  or  H,^^  ~  uH^. 
In  fetch  dependent  case,  an  empirical  equation  for  Tm  is  proposed  here  as 

§^^d,['-^•t^nh{d,y)\'^  (4.16) 

where  di,  ^2,  and  dz  are  constants.  For  large  F/F,,  Eq.(4.16)  approciches  to 

W^  =  M^J^r\    (orF  »   F,  (4.17) 

In  duration  dependent  case,  Tm  is  proposed  here  cis 

^  =  rf.[^-tanh(.isf)]^«  (4.18) 
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where  d^,  d^  and  de  are  constants.  For  large  t/t,,  Eq.(4.18)  approaches  to 

'-^-dA'^,     iort»U  (4.19) 

Thus,  based  on  either  Eq.(4.17)  or  (4.19),  the  upper  limit  of  Tm  corresponding  to 
Uio  can  be  approximated  by  the  following  expression: 

^^  =  dr{U,oY'  (4.20) 

where  Tm,oo  indicates  the  modal  period  at  the  fully  developed  sea.  If  this  fully 
developed  condition  is  approximated  by  a  upper  bound  of  Tm  —  Uio  curves,  drawn  as 
a  solid  line  in  Fig.  4.15,  we  find  dj  w  43.0  (m/sec)^/^  and  ds  «  -2/3,  or  Tm,oo  ~  U^^. 
In  order  to  exajnine  the  relationship  between  the  resonant  parameter,  2-kUiqI gTm 
(=  fm/fo),  and  the  energy  scale  parameter,  Hg/gT^  (=  f^HJg),  when  the  wind 
generated  waves  develop  into  a  saturation  condition,  their  computed  values  bcLsed 
on  the  UCWP  model  are  plotted  against  each  other  in  Fig.  4.16.  It  can  be  seen  that 
the  fetch  limited  case  and  the  duration  limited  case  behave  in  a  similar  manner  and 
both  approach  to  the  same  limiting  curve.  The  long-dash  lines  shown  in  Fig.  4.16 
indicate  the  fully  developed  condition  solved  from  Eqs.(4.12)  and  (4.17)  by  elimi- 
nating the  common  term  gF,/U^Q,  or  from  Eqs.(4.14)  and  (4.19)  by  eliminating  the 
common  term  gtg/Uio: 


2-kU,.        H, 


Also  shown  in  Fig.  4.16  is  the  saturation  condition,  drawn  in  short-dcish  lines, 
obtained  from  the  classical  SMB  method  which  utilizes  [/„  =  Uiq  in  Eqs.(4.1)  and 
(4.6): 

It  is  seen  that  the  UCWP  model  yields  greater  energy  scale  parameter  than  the 
SMB  model  when  the  resonant  parameter  is  small  (approximately  less  than  unity), 
and  it  yields  less  energy  scale  parameter  than  the  SMB  model  when  the  resoneint 
parameter  is  large  (approximately  greater  than  unity). 
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Figure  4.16:  Plots  of  27rf/,o/yT^  versus  lOHlJgTl  'n  (a)  fetch-limited  case  and  in 
(b)  duration-limited  czise. 
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The  resonajit  parameters  ajid  the  energy  scale  parameters  baised  on  measured 

buoy  data  are  then  used  to  show  the  contour  lines  of  Hg  in  Fig.  4.17.  Also  shown  are 
the  saturation  lines  given  by  the  UCWP  model  eind  the  SMB  model.  These  limits 
are  good  indicators  of  the  maximal  wave  height  imder  optimum  wind  conditions.  It 
is  seen  that  the  UCWP  model  shows  better  agreement  with  the  measured  data  thaji 
the  classical  SMB  method  for  high  wave  condition.  For  instance,  the  long-dash  line 
which  represents  the  saturation  relation  obtained  from  the  UCWP  model  passes 
through  the  centers  of  the  maximum  contoxir  of  7-meter  for  the  1984  data  and  the 
one  of  8-meter  for  the  1985  data  whereas  the  short-dash  line  which  represents  the 
classical  SMB  model  crosses  the  maximum  contour  only  for  the  1984  data. 

Waves  above  the  saturation  condition  in  Fig.  4.17  are  in  the  growing  stage  and 
they  are  influenced  by  the  limited  fetch  and  duration.  Waves  below  the  long-dash 
line  are  in  the  decaying  stage  when  experiencing  the  dissipation  under  the  relatively 
low  winds. 

4.5     Comparison  of  Various  Wave  Energy  Fliix  Elements 

For  wave  generation  in  the  deep  water,  the  turbulent  dissipation  has  been  men- 
tioned to  be  a  dominant  factor  in  the  UCWP  model.  Its  contribution  to  the  total 
energy  flux  is  to  be  compared  with  the  advection  and  wind  input  terms  in  this  sec- 
tion. The  comparisons  are  shown  for  three  cases;  the  first  two  consider  the  growing 
and  decaying  stages,  and  the  last  one,  the  small  wind  and  wave  condition.  The  first 
case.  Case  1,  is  during  a  growing  period  under  increasing  wind  condition  at  1:00  a.m. 
September  8,  1984,  when  the  winds  increasing  from  26  knots  to  31  knots.  A  net 
gain  of  wave  energy  is  observed  in  the  case.  The  second  case,  Case  2,  is  a  decaying 
case  at  5:00  a.m.  September  8,  1984,  while  onshore  winds  decreased  from  30  knots 
to  28  knots  in  two  hours.  A  net  loss  of  energy  is  noticed  in  the  case.  The  last  case. 
Case  3,  shown  in  Fig.  4.20  is  a  loczd  wind  wave  growth  case  at  23:00  p.m.  October 
13,  1984,  when  winds  remain  a  constant  speed  at  12  knots  for  the  last  twelve  hours. 
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Figure  4.17:    Plots  of  27rf/,o/^r„  versus  lO^HjgTl  based  on  the  NOAA  buoy 
#41006  1984-1985  data. 
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Table  4.3:  Summary  of  energy  fliix  budget  at  the  buoy  #41006  location. 


mechanism 

lO^AE/At     (mVsec) 

Case  1 

Case  2 

Case  3 

advection 

-2.93 

-4.80 

-0.09 

Gin,i  wind  input  term 

13.97 

3.01 

0.64 

Gin,2  wind  input  term 

0.63 

0.64 

0.08 

turbulent  viscous  loss 

-3.91 

-5.15 

-0.65 

net  (temporal  change) 

7.76 

-6.30 

-0.02 

Table  4.3  shows  the  energy  flux  budget  for  the  three  cases  at  the  buoy  #41006  loca- 
tion. In  all  three  cases,  the  turbulent  dissipation  is  important.  The  local  wind  input 
term,  G,n,2  is  overall  relatively  small  in  all  three  cases.  In  Case  3,  where  waves  are 
nearly  in  stationary  condition  under  constant  wind  speed,  temporal  chcinge  is  found 
to  be  very  small  when  the  G.^,!  term  is  approximately  balanced  by  the  turbulent 
dissipation  term.  Plots  of  rate  of  change  of  individual  deepwater  effect  sources  over 
the  frequencies  for  all  three  cases  are  given  in  Appendix  D. 

It  is  found  here  that  the  bottom  effects  in  shallow  water  are  quite  pronounced 
for  the  Florida  Coast.  Three  cases  showing  the  contributions  of  various  elements  to 
the  energy  flux  at  a  shallow  water  location  are  presented  here.  In  each  case  the  rates 
of  change  of  spectrum  due  to  the  temporal  variation,  the  advection,  which  includes 
shoaling  and  refraction  effects,  the  atmospheric  input,  the  turbulent  dissipation, 
and  the  bottom  frictional  loss,  are  evaluated  from  the  shallow  water  UCWP  model. 
Thus,  the  relative  importance  of  individual  spectral  components  can  be  noticed 
directly  from  these  cases. 

The  first  case.  Case  4,  is  during  a  growing  period  under  increasing  wind  condi- 
tion at  5:00  a.m.  September  8,  1984,  when  the  onshore  winds  increasing  from  27 
knots  to  30  knots  in  two  hours.  A  net  gain  of  wave  energy  is  observed.  The  second 
case,  Case  5,  is  a  decaying  case  at  1:00  a.m.  September  10,  1984,  while  onshore 
winds  decreased  from  34  knots  to  29  knots  in  two  hours.    A  net  loss  of  energy  is 
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Table  4.4:  Siimmajy  of  energy  fliix  budget  at  the  Marineland  station. 


mechanism 

lO^AE/At    (mVsec) 

Case  4 

Case  5 

Case  6 

advection,  shoaling, 
ajid  refraction 

1.79 

1.75 

1.68 

atmospheric  input 

13.72 

3.13 

3.72 

turbulent  viscous  loss 

-0.77 

-0.68 

-0.45 

bottom  frictionaJ  loss 

-9.15 

-11.07 

-4.41 

net  (temporal  change) 

5.59 

-6.87 

0.54 

noticed  in  the  case.  The  third  case,  Case  6,  occtirred  at  3:00  a.m.  September  10, 
1984,  when  onshore  winds  were  dropping  from  29  knots  to  26  knots  in  two  hours. 
In  this  last  case,  a  net  loss  of  energy  is  seen  at  low  frequencies  and  a  net  gain  of 
energy  is  found  at  the  high  frequency  range,  which  is  due  to  the  entering  of  high 
frequency  waves  from  the  neighboured  area.  Table  4.4  shows  the  energy  flux  budget 
for  the  three  cases  at  the  Marineland  station.  In  the  all  three  cases,  the  bottom 
friction  and  advection  effects  dominate  the  other  processes.  More  than  65  percent 
of  the  wind  input  is  being  dissipated  by  bottom  friction.  Temporal  changes,  often 
neglected  in  shallow  water  wave  transformation  computations,  are  found  to  be  also 
significant.  The  dissipation  of  wave  energy  due  to  turbulent  viscosity  is  seen  to 
be  less  important  in  these  cases.  The  plots  of  rate  of  change  of  individual  shallow 
water  effect  sources  over  the  frequencies  for  the  above  three  cases  are  also  given  in 
Appendix  D. 

4.6     Influence  of  c/  in  Shallow  Water  Wave  Transformation  Process 

An  important  factor  that  influences  the  shallow  water  wave  transformation  is 
the  bottom  friction.  To  estimate  it  according  to  Eq.(2.50)  in  a  wave  model,  bottom 
friction  coefficient  Cf  needs  to  be  first  determined.  The  magitude  of  Cf  has  been 
suggested  to  range  from  0.005  to  0.05  at  various  locations  along  the  Florida  coast. 
Five  values  of  the  friction  coefficient,  Cf  =0.0,  0.005,  0.01,  0.02,  and  0.04,  are  selected 
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here  to  test  their  influences  on  the  bottom  frictional  loss.  Figure  4.18  shows  the 
comparsions  of  wave  heights  and  modal  periods  hindcasted  in  September  1984  with 
respect  to  the  five  Cj  values.  Both  the  time  series  of  the  hindcasted  wave  height 
and  modal  periods  plotted  in  Fig.  4.18  were  smoothed  by  three  points  running 
average  over  a  six-hour  interval  to  provide  better  visual  comparison.  As  expected, 
waves  hindcasted  from  the  UCWP  model  based  on  C/=0.0  have  the  largest  heights 
and  modal  periods  than  the  other  four  cj  values.  That  is,  when  the  dissipation 
mechanism  due  to  bottom  friction  is  absent  in  the  model,  long  waves  can  survive 
and  result  more  wave  energy  in  the  low  frequency  region.  On  the  other  hands,  waves 
hindcasted  based  on  C/  =  0.04,  which  is  the  largest  among  the  five  tested  values,  are 
seen  to  generate  the  least  heights  ajid  modal  periods  compared  with  the  ones  based 
on  the  other  four  Cj  values.  This  is  because  tremendous  energy  has  been  dissipated 
in  the  low  frequency  regions  due  to  the  utilization  of  strong  bottom  friction  in  the 
model.  The  wave  heights  and  modal  periods  computed  from  the  model  based  on 
the  c/  equal  to  0.005,  0.01,  and  0.02,  respectively,  are  seen  to  be  monotonically 
decreasing  between  the  ones  with  c/=0.0  (null  bottom  friction  condition)  and  the 
ones  with  c/=0.04  (strong  bottom  friction  condition).  By  using  different  Cy  values 
in  the  UCWP  model,  the  wave  heights  are  seen  to  change  siginificantly  for  large 
wave  condition  but  insignificantly  for  small  wave  condition.  However,  the  modal 
periods  may  change  dramatically  for  small  but  long  waves  since  the  bottom  frictional 
dissipation  is  more  active  in  low  frequency  components  than  in  high  frequency 
components.  The  change  of  modal  period  for  high  waves  according  to  various  c/ 
values  is  relatively  small  but  detectable. 

Also  included  in  Fig.  4.18  are  the  meaisured  wave  heights  and  modal  periods 
which  allow  a  statistical  analysis  between  the  computed  and  measured  data.  The 
statistic  parameters  chosen  here  for  the  analysis  are  the  correlation  coefficient  and 
a  root-mean-square  error  normalized  by  the  root-mean-square  value  of  measured 
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Figure  4.18:  Comparisons  of  wave  heights  and  modal  periods  computed  in  Septem- 
ber 1984  for  c/=0.0,  0.005,  0.01,  0.02,  and  0.04. 
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Table  4.5:    Comparison  of  <5„.,,„..„,  <5^„,,r^,„,  £«.,,«.,„,  and  e 
0.005,  0.01,  0.02,  and  0.04. 


Tm.,e  Tm,m 


for  c/=0.0, 


<'/ 

6' 

H..cH,,m 

6' 

Tm,c  Tfn.m. 

H»,c  H»,m 

Tm,c  Tm,m. 

0.0 

0.93 

0.30 

0.43 

0.52 

0.005 

0.94 

0.32 

0.30 

0.51 

0.01 

0.94 

0.31 

0.25 

0.51 

0.02 

0.95 

0.38 

0.23 

0.37 

0.04 

0.92 

0.30 

0.26 

0.51 

*Computations  are  based  on  the  Sep.  1984  data  at  Marineland  Station 


data.  The  latter  can  be  expressed  as: 


t=l 


1  t=l 


(4.23) 


where  X^  is  a  computed  X,  Xm  is  a  measured  X,  and  N  is  the  sample  size.  The 
higher  correlation  and  the  lower  normalized  root-mean-square  error  indicate  better 
agreement  between  the  computed  and  measured  data.  Table  4.5  shows  the  compar- 
ion  of  these  statistical  parameters  between  computed  and  measured  wave  heights 
and  modal  period  for  c/=0.0,  0.005,  0.01,  0.02,  and  0.04.  It  is  noticed  the  cor- 
relations between  the  computed  and  measured  wave  heights  based  upon  the  five 
different  c/  values  are  all  high  and  they  do  not  differ  significantly  from  one  to  the 
others  although  the  one  based  on  c/=0.02  gives  the  hightest  correlation.  The  cor- 
relations between  the  computed  and  measured  modal  periods  are  found  between 
0.3  and  0.38  based  on  the  five  c/  values.  The  one  based  on  c/=0.02  again  gives 
the  highest  correlation.  The  comparison  of  normalized  root-mean-square  error  for 
different  Cj  values  shows  that  errors  calculated  from  wave  height  and  modal  period, 
respectively,  based  on  c/=0.02  Jire  smaller  than  the  others.  Therefore,  the  utiliza- 
tion of  c/=0.02  in  computing  bottom  friction  is  good  for  the  wave  hindcast  along 
the  Florida  coast. 

Figure  4.19  (a),   (b),   (c),   (d),  and  (e),  respectively,  give  examples  of  how 
the  bottom  frictlonal  dissipation  contributes  to  the  total  energy  flux  for  C/=0.0, 
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Figure  4.19:  (a)  Plot  of  individual  shallow  water  effect  sources  for  c/=0.0 
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Figure  4.19:  (b)  Plot  of  individual  shallow  water  effect  sources  for  c/=0.005 
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Figiire  4.19:  (c)  Plot  of  individual  shallow  water  effect  sources  for  c/=0.01 
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Figure  4.19:  (d)  Plot  of  individual  shallow  water  effect  sources  for  c/=0.02 
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Figure  4.19:  (e)  Plot  of  individual  shallow  water  effect  sources  for  C/=0.04 
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Table  4.6:  Summary  of  energy  flux  budget  for  cases  with  c/=0.0,  0.005  0.01,  0.02, 
and  0.04,  respectively. 


mechanism 

10^A£;/Ai     (mVsec) 

Cf  =  0 

Cf  =  0.005 

Cf  =  0.01 

Cf  =  0.02 

Cf  =  0.04 

advection,  shoaling, 
and  refraction 

-1.59 

0.00 

0.89 

1.79 

2.02 

atmospheric  input 

6.86 

10.27 

11.77 

13.72 

17.63 

tiirbulent  viscous  loss 

-1.25 

-1.06 

-0.90 

-0.77 

-0.70 

bottom  frictional  loss 

0.00 

-4.66 

-6.77 

-9.15 

-13.13 

net  (temporal  change) 

4.02 

4.55 

4.99 

5.59 

5.82 

0.005,  0.01,  0.02,  and  0.04  at  5:00  a.m.  September  8,  1984.  The  bottom  frictional 
dissipation  increases  from^  zero  at  no  bottom  friction  situation  in  Fig.  4.19(a)  to 
the  greater  frictional  loss  cases  in  Figs.  4.19(b),  (c),  (d),  and  (e).  Table  4.6  shows 
the  summary  of  energy  flux  budget  for  the  five  caises  indicated  in  Fig.  4.19.  In 
stronger  bottom  friction  caise  which  correponds  to  higher  Cf  value,  more  energy 
flux  components  associated  with  wind  input  are  seen  to  shift  to  the  high  frequency 
region.  This  is  because  in  the  stronger  friction  case  more  energies  are  dissipated  in 
the  low  frequency  region  than  in  the  high  frequency  region. 


CHAPTER  5 
CONCLUSIONS 


5.1     Summaxy 

This  study  examines  wave  models  in  hindcasting  waves  in  an  eixea  of  the  Atlantic 
Ocean  covering  the  southecistem  coast  of  the  United  States.  For  wave  models,  the 
input  required  is  the  wind  information  over  the  sea  surface.  The  numerical  methods 
employed  in  wave  models  cire  based  upon  the  energy  balance  equation  which  states 
that  the  total  difference  of  wave  energy  is  equal  to  the  summation  of  the  energy 
fluxes  due  to  wave  growth,  dissipation,  and  wave-wave  interaction.  Wave  models  in 
general  consist  of  a  deepwater  model  ajid  a  shallow  water  model.  The  shallow  water 
model  includes  the  influences  from  the  shallow  water  effects,  such  as  refraction, 
shoaling,  percolation,  bottom  friction,  etc.,  which  are  excluded  in  the  deepwater 
model. 

In  deep  water,  the  growth  of  waves,  as  a  result  of  atmospheric  input,  is  a  highly 
nonlinear  process.  The  growth  rate  associated  with  this  mechanism  is  usually  es- 
timated in  deepwater  model  from  two  contributions:  (i)  the  growth  of  prevailing 
waves,  which  is  tziken  to  be  proportional  to  the  third  power  of  am  energy  scale  pa- 
rameter, a,  and  (ii)  the  growth  of  local  wind  waves,  which  is  taken  to  be  linearly 
proportional  to  a.  The  dissipation  of  waves  is  mainly  from  wave  breakings  which 
generates  turbulence  induced  in  wave  motion.  This  dissipation  mechanism  is  gen- 
erally considered  in  the  model  by  maintaining  an  f~^  law  on  the  down-slope  face  of 
spectrum.  The  advection  is  often  transacted  in  the  model  by  using  a  free  propaga- 
tion algorithm.  Thus,  wave  energies  between  directional  and  frequency  bands  are 
not  coupled. 
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In  shallow  water,  the  refraction,  shoaling,  percolation,  and  bottom  friction  be- 
come important  in  the  wave  energy  transport  process.  Wave  paths  are  no  longer 
independent  on  water  depth.  Wave  rays  axe  traditionally  determined  in  the  shal- 
low wave  model  from  the  equation  of  nondivergence  of  wave  nimiber  vectors  in  the 
horizontal  plane.  Shoaling  ajid  refraction  are  taJien  into  account  in  the  model  by 
redistributing  wave  energies  ajnong  the  neighbouring  wave  rays  according  to  the 
wave  balance  equation.  There  is  no  coupling  of  wave  energies  among  individual 
frequency  bands.  The  percolation  and  bottom  friction  effects  are  considered  as  sink 
terms  in  the  model.  The  spectrzil  growth  rate  due  to  wind  input  in  the  model  com- 
monly follows  a  combined  Phillips'  ajid  Miles'  growth  mechanism,  which  says  that 
the  growth  rate  is  a  linear  function  of  wave  spectnma. 

Two  previouly  developed  deepwater  wave  hindcasting  models,  the  discrete  model 
by  Wave  Information  Study  (WIS)  group  (Resio  and  Tracy,  1983)  and  the  parzimet- 
ric  model  by  Hasselmann  et  al.(l973),  are  reviewed.  It  is  demonstrated  in  the 
present  study  that  in  fetch-limited  cases  the  two  models  generally  consistent  in  pre- 
dicting wave  heights  and  modal  periods  for  small  fetch  (less  than  400  km),  and 
the  parametric  model  predicts  much  higher  wave  heights  and  model  periods  than 
the  WIS  model  for  large  feth  (greater  than  400  km).  In  duration-limited  cases  the 
two  models  are  consistent  in  predicting  the  significant  wave  heights  only  when  the 
duration  is  less  than  20  hours.  The  parametric  model  predicts  again  much  higher 
wave  heights  and  modal  periods  than  the  WIS  model. 

The  new  deepwater  model  developed  here  is  of  discrete  type  and  is  based  on 
WIS  model  with  modified  growth  and  dissipation  terms.  The  wave  growth  rate  due 
to  atmospheric  input  is  again  estimated  in  the  model  for  the  growths  of  prevailing 
waves  and  local  wind  waves.  The  growth  rate  of  prevailing  waves  are  btised  upon  the 
modified  JONSWAP  spectrum  defined  by  Donelan  et  al.(l985)  and  it  is  proportional 
to  a'  which  reflects  the  nonlinear  nature  in  accounting  for  the  wave-wave  interaction. 
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The  growth  of  local  wind  waves  is,  however,  taken  to  be  one  of  the  three  cases:  (a)  if 
the  shift  of  wind  direction  from  the  previous  time  step  is  quite  large,  new  fresh  wind 
waves  will  be  generated,  (b)  if  the  shift  of  wind  direction  from  the  previous  time  step 
is  small  but  not  null,  wind  waves  generated  from  the  previoiis  time  step  will  continue 
to  grow  in  the  direction  of  the  wind,  and  (c)  if  the  shift  of  wind  direction  is  small 
from  the  previous  time  step  and  the  wind  direction  is  the  same  as  the  prevailing 
wave  direction,  this  source  term  will  vanish.  In  the  first  case  the  generation  of  new 
waves  is  based  on  the  spectral  form  of  P-M  spectrum  and  is  proportional  to  a  which, 
therefore,  reflects  a  linear  growth  process.  In  the  second  case  the  growth  of  wind 
waves  is  treated  as  a  minor  nonlinear  process  which  is  proportional  to  a^.  In  the 
last  case  the  growth  of  wind  waves  follows  the  same  mechanism  as  the  growth  of 
prevailing  waves.  The  turbulent  dissipation  is  estimated  in  the  model  according  to 
the  algorithm  proposed  by  Hasselmann  (1974).  The  dissipation  rate  is  expressed  as 
the  product  of  the  spectrum,  a  wave  frequency,  and  a  proportional  constant.  The 
proportional  constant  is  calibrated  by  means  of  field  ocean  wave  spectral  data. 

The  test  of  the  new  deepwater  model  is  performed  in  both  hypothetical  and  real 
cases.  In  the  hypothetical  cases,  the  new  model  is  tested  against  the  WIS  model, 
the  SMB  model  and  the  parametric  model.  In  fecth-limited  cases  the  significant 
wave  heights  obtained  from  the  model  agree  well  with  those  obtained  from  the  WIS 
model  at  small  winds  (20  and  30  knots),  and  they  agree  well  with  those  obtained 
from  the  SMB  model  at  moderate  wind  speeds  (40  knots);  they  are  higher  than 
those  obtained  from  the  other  three  models  at  high  wind  speeds  (50  knots)  when 
fetch  is  less  than  400  km.  The  modal  periods  obtained  from  the  new  model  are 
consistent  with  those  obtained  from  the  SMB  model  at  low  wind  speeds  (20  and  30 
knots).  At  higher  wind  speed,  the  modal  periods  obtained  from  the  new  model  are 
much  higher  than  the  other  three  models. 
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In  the  hypothetical  duration-limited  cases  the  significant  wave  heights  obtained 
from  the  new  model  agree  well  to  those  obtained  from  the  SMB  model  when  du- 
ration is  less  than  20  hrs.  For  duration  greater  than  20  hrs,  the  significant  heights 
computed  from  the  new  model  are  not  much  different  from  those  computed  from 
the  WIS  model  except  at  40  knots  winds  where  wave  heights  obtained  from  the 
WIS  model  has  a  non-realistic  damping  pattern.  The  wave  heights  obtained  from 
the  new  model  are  much  smaller  than  those  obtained  from  the  SMB  and  paraunetric 
models  for  duration  greater  than  20  hrs.  The  model  periods  computed  from  the 
new  model  are  not  much  different  from  those  computed  from  the  SMB  model  at 
small  and  moderate  wind  speeds  (20,  30,  and  40  knots);  they  are,  however,  much 
smaller  than  those  obtained  from  the  parametric  model  when  diiration  is  greater 
than  20  hrs. 

In  the  real  cases,  the  new  deepwater  model  is  compared  with  measured  waves. 
The  wave  heights  computed  from  the  model  are  highly  correlated  to  those  measured 
by  NOAA  maintained  buoys.  The  modal  periods  computed  from  the  model  are  not 
well  correlated  to  the  measured  ones.  However,  the  computed  modal  periods  are 
found  to  be  more  in  accord  with  the  measured  data  for  large  waves.  Swells  not  be 
properly  traced  in  the  model. 

The  shallow  water  model  developed  in  this  study  follows  Chen  and  Wang's 
(1983)  model  with  modified  generation  term.  The  refraction  and  shoaling  are  eval- 
uated with  the  advection  term  in  the  energy  balance  equation.  Percolation  and 
bottom  friction  are  included  in  the  model  by  using  the  formulas  derived  by  Shemdin 
et  al.(1978)  and  Collins  (1972),  respectively.  The  turbulent  dissipation  is  estimated 
in  the  model  the  same  as  in  the  deepwater  model.  The  spectral  growth  rate  due  to 
wind  input  is  estimated  by  using  the  same  functional  form  for  local  wave  genera- 
tions in  the  deepwater  case;  it  depends  not  only  on  the  strength  of  the  current  wind 
strength  but  also  on  the  current  wave  condition. 
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5.2     Conclusions 

A  coupled  discrete  spectral  model,  named  here  as  the  Unified  Coastal  Wave 
Prediction  (UCWP)  model,  has  been  developed  using  surface  pressure  information 
as  input  to  predict  wind  waves  in  both  deep  and  shallow  water  regions.  The  model 
has  been  applied  to  hindcasting  wind  waves  along  the  east  coast  of  Florida  and  it 
performed  well  for  both  low-  and  high-pressure  weather  systems.  The  significant 
wave  heights  ajid  directions  of  wave  propagation  are  highly  correlated  to  the  mea- 
sured data.  The  results  of  modal  period  prediction  are  good  only  for  high  wind 
waves  but  poor  for  low  wave  conditions. 

The  fimdamental  processes  employed  in  the  UCWP  model  are  similar  to  those 
used  by  other  widely-accepted  discrete  and  paraxaetric  models.  New  source  and 
sink  terms  are  proposed  on  the  basis  of  field  evidence. 

Several  conclusions  relating  to  this  research  are  summarized  in  the  following 
paragraphs. 

(l)  Plots  of  spectral  peak  energy  versus  spectral  peak  frequency  in  the  wave 
growth  cases  show  that  the  data  are  very  scattered  in  the  equilibrium  range  indi- 
cated by  the  P-M  spectrum.  Although  the  data  scattering  could  be  exaggerated 
by  the  presence  of  swell,  which  generally  appears  as  long  waves,  there  is  no  indica- 
tion that  a  simple  ciirve,  such  as  the  P-M  equilibrium  spectrum,  can  represent  the 
relationship.  The  scattering  of  the  data  suggests  active  nonlinear  wave-wave  inter- 
actions. That  is,  energy  transfer  from  higher  frequency  components  to  lower  fre- 
quency components  is  a  nonlinear  process  which  is  dependent  upon  multi-variables 
including  wind  speed,  wave  phase  speed,  peak  energy  frequency,  etc.  The  nonlinear 
energy  transfer  can  be  even  stronger  as  indicated  by  the  plots  of  nonlinear  trans- 
fer coefficient  Uni  versus  /o//^Af  which  show  that  a^  is  more  likely  a  function  of 
/o//^Af  than  a  constant  value  suggested  by  Hasselmann  (1976)  and  used  in  the 
WIS  deepwater  model. 


120 

(2)  In  the  growth  of  prevailing  wave  components,  energy  transfers  due  to  atmo- 
spheric input  and  wave-wave  interactions  are  both  highly  nonlinear  and  of  similar 
nature.  Thus,  they  can  be  treated  as  one  growth  term.  It  also  implies  that  the 
energy  transfer  process  is  mainly  governed  by  the  components  receiving  the  en- 
ergy irrespective  of  energy  source  forms.  The  scale  paramter,  a,  which  controls  the 
nonlinear  growth  rate  is  approximated  by  a  power  function  of  fmal  fo-  A  power  of 
7/6  is  selected  based  on  data,  instead  of  the  value  of  2/3  as  suggested  by  Hzissel- 
mann  et  al.(l976).  This  results  in  a  much  stronger  growth  under  severe  weather. 

(3)  The  'overshoot'  phenomenon,  which  indicates  waves  with  slightly  higher 
height  than  the  saturation  height  occur  before  reaching  a  fully-develped  condition, 
is  observed  in  the  test  of  the  UCWP  deepwater  model  in  the  duration-limited  case. 
The  phenomenon  is  in  general  quite  weak  as  revealed  by  the  data  but  essential  in 
building  enough  energy  to  reach  the  later  saturation  condition.  It  is  not  found  in 
fetch-limited  case,  and  it  is  neither  observed  in  the  WIS  model  nor  in  the  SMB  and 
parametric  models. 

(4)  Based  upon  wave  data  obtained  from  the  UCWP  deepwater  model,  both 
significant  wave  height  ajid  modal  period  can  be  approximated  by  power  functions 
of  wind  speed  at  infinite  fetch  and  duration.  The  significant  height  is  found  to 
be  proportional  to  wind  speed  with  the  power  equal  to  5/3  and  the  modal  period 
is  found  to  be  proportional  to  wind  speed  with  the  power  equal  to  1/3.  These 
two  empirical  equations  were  then  combined  to  yield  a  fully-developed  condition 
and  compared  with  the  field  wave  data.  The  comparison  shows  that  the  combined 
empirical  relation  agrees  well  to  these  measured  data  at  high  wave  condition.  The 
fully-developed  condition  given  by  the  classical  SMB  model  is  also  plotted  against 
the  same  field  data  set.  It  is  shown  that  the  fully-developed  relation  between  wind 
and  wave  obtained  from  the  classical  SMB  model  is  less  consistent  with  the  measured 
data  than  the  one  obtained  from  the  UCWP  deepwater  model. 
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(5)  Dissipation  of  wave  energy  due  to  tiirbulence  included  in  the  UCWP  model 

follows  Hasselmajin's  formula  which  is  proportional  to  the  product  of  a  wave  spec- 
trum and  wave  frequency.  Its  contribution  to  the  total  energy  flux  is  found  to  be 
important  compared  with  the  advection  and  wind  input  terms  in  the  deepwater 
spectral  tranformation  process.  However,  it  is  generally  of  a  small  order  of  magni- 
tude compared  with  the  advection,  shoaling,  refraction,  and  bottom  friction  terms 
in  the  shallow  water  spectral  transformation  process. 

(6)  In  shallow  water,  bottom  friction  always  plays  aji  important  role.  It  con- 
sumes much  of  the  wind  energy  input  for  the  cases  tested  along  the  Florida  coast. 
Long  swell  components  are  severely  curtailed  when  propagating  from  deep  to  shal- 
low water.  The  waves  in  shallow  water  could  be  highly  nonstationary  and  this  is 
particularly  the  case  under  high  winds.  Hence,  the  validity  of  ray  method  which 
can  only  be  applied  to  stationary  cases  is  doubtful  for  high  wind  waves. 

(7)  Test  of  response  of  the  UCWP  shallow  water  model  to  diff"erent  values  of 
bottom  friction  coefficient,  c/,  shows  some  interesting  results.  By  using  diff"erent 
c/  values  in  the  UCWP  model,  the  wave  heights  caji  change  siginificantly  for  large 
wave  condition  but  insignificantly  for  small  wave  condition.  The  modal  periods  may 
change  dramatically  for  small  but  long  waves  since  the  bottom  frictional  dissipation 
is  more  active  in  low  frequency  components  than  in  high  frequency  components.  The 
change  of  modal  period  for  high  wave  condition  according  to  various  Cf  values  is  in 
general  small.  The  statistical  analysis  between  the  computed  and  measured  wave 
heights  based  upon  diff"erent  Cj  values  shows  that  correlation  coefficients  are  overall 
high  and  they  do  not  diff"er  significantly  from  one  to  the  others  although  the  one 
based  on  c/=0.02  gives  the  hightest  value.  The  correlations  between  the  computed 
cind  measured  modal  periods  are  found  between  0.3  and  0.38  for  all  the  test  cases;  the 
one  based  on  c/=0.02  again  gives  the  highest  value.  The  comparison  of  normalized 
root-mean-square  error  for  diff"erent  cj  values  shows  that  errors  calculated  from 
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wave  height  and  modal  period,  respectively,  based  upon  c/=0.02  are  smaller  thaji 
the  others.    Therefore,  Cf— 0.02  is  chosen  in  the  UCWP  shallow  water  model  for 
wave  hindcast  along  the  Florida  coast. 

5.3     Future  Studies 

Further  studies  related  to  this  work  are  suggested  as  follows: 

(1)  Experiments  producing  deepwater  waves  in  the  ideal  wind  condition  need 
to  be  conducted  in  a  large  scale  physical  model  in  order  to  establish  a  dependable 
data  beise  for  the  use  in  calibration  of  numerical  wave  models. 

(2)  Wave-current  interaction  effect,  which  is  definitely  important  in  the  presence 
of  strong  current  in  the  wave  field,  should  be  included  in  the  wave  model. 

(3)  Performance  of  a  deepwater  wave  model  in  the  entire  ocean  basin  is  required 
so  as  to  include  the  complete  swell  system  and,  therefore,  to  increase  the  ability  in 
predicting  the  peak  energy  frequency  as  well  as  the  significant  wave  height. 


APPENDIX  A 
SUMMARY  OF  PARAMETERIZED  PROGNOSTIC  EQUATIONS 


A.l     Prognostic  Equations 

dbi  dbi 

-^  +  A;— ■  =  G,-,  {,j  =  1,2,3,4,5  (A.l) 

where  61  =  fm,  62  =  a,  63  =  7,  64  =  (7a,  65  =  ctj,  i.  represents  the  horizontal  axis 
along  which  waves  propagate,  I?,y  denotes  the  matrix  of  generalized  propagation 
velocity,  ajid  G,  denotes  the  vector  of  projected  soiirce  function. 

A. 2     Elements  of  Z),y  and  G, 


D,,  =  v{l  +  5K),  D,2  =  -v~K,  Di3  =  -v^K, 

a  7 

r>2i  =  t;-^(0.263  +  3.61K),       D22  =  u(0.605  -  0.722is:),    D23  =  -V-0.722K, 

fm  1 

D2i  =  -V ^{0.263-21  AK),    D32  =  v-{OA  -  A.28K),       D33  =  v(l  -  4.27SK), 
fm  a. 

D41  =  -^{5K{aa  -  1)  +  cTa  +  7^[8.2436<7„  +  0.1705  -  /f (164.87(7^  +  13.877)]}, 
fm  In  7 

Di2  =  -{K(<7a  -  1)  +  T^(-A"(2.7755  +  32.974<7J  +  1.6487c7a  +  0.25644]}, 
a  ln7 

As  =  -{K{oa  -  1)  +  r^[<7a  -  i^(32.974(7„  +  2.7755)]}, 
7  In  7 

Ai  =  --^{Si^K  +  1)  +  (76  -  -^[8.2436a6  -  0.1705  -  K{ie4.87(j^  -  13.877)]}, 
fm  In  7 

A2  =  -{i<'(<76  +  1)  -  p-(K(2.7755  -  32.974afc)  +  1.6487cr6  -  0.25644]}, 

^3  =  -{K{a,>  +  1)  -  T^kj  -  A'(32.974a4  -  2.7755)]} 
7  ln7  ^ 

(A.2) 
where  v  =  fir/4  tt/^,  A"  =  a^/{2Qa'^  +  ln7),  and  ct  =  (a^  +  <76)/2. 
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Gi  = 


-0.586  aV^^ 


0, 


,       if  ^  >  0.13, 


otherwise, 


G2  =  (0.005022  t^^  -5a^)afm, 
G3  =  16(70  -  1)  ct^fm, 
G4  =  [25.5(0.07x  -  aa)  -  0.5(^6  -  0.09x)]  aVm, 
G5  =  [25.5(0.09x  -  aj)  -  0.5(a<,  -  0.07x)]  a'/m, 


(A.3) 

(A.4) 
(A.5) 
(A.6) 

(A.7) 


where 


and 


e  = 


/mf/i 


10 


^       ^  +  0.7^ 


3.3,  if  ^  >  0.13, 

£-0.13 
'  if  0.16  >   £  >  0.13, 


1  +  2.3- 
1, 


0.03 


if  0.13  >   ^ 


(A.8) 


(A.9) 


APPENDIX  B 
PROCEDURE  IN  SOLVING  PARAMETERIZED  G,„,i(/) 


a..AM  =  j^J  ^M-(^/]  (B.i) 


where 


(Ta,         for  /n    <    frri, 

<^ab=  \  n=  1,2,  •••,n,„,---,N 

<^6,         for  /n    >  /^, 


Urn  —  index  of  frequency  at  maximum  E{fn), 
fm  =  frequency  at  maximum  E{fn), 
N  =  total  number  of  frequency  bands 

The  parameters  are  determined  by  minimizing  the  root-mean-square  error 


ElG'.n.lC/nl-G'C/JlVEGC/J  (B.2) 

n=l 


where 

G(/„)  =  AEiU)/At 

is  obtained  based  on  measured  wave  spectra.  Taking  the  derivatives  of  (rma  with 
respect  to  a,  7,  tr,,,  c^  and  equating  them  to  zero,  four  equations  can  be  obtained 
and  expressed  as  follows: 

-'  Y.-~^-M-'^{fY\  =  E  Wn)^exp(-(^n  (B.3) 
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-'  E  ^^TTw-^^P'-'fr'*!  =  E£{/n)^exp|-(^)*l  (B.4) 

n=l  \^^)    JnJm  Jn  „=!  J n  •'n 

nm      „2o      20.  / 

n=l  Ja  Jn 


=   E  £(«^(/. -/■»)' exp|-(^)<]  (B.6) 

n=nn»  Jn  Jn 

where  /?a  and  /3j  aire  corresponding  to  <7a  and  aj,  respectively. 

In  order  to  solve  these  non-linear  equations,  a  numerical  iteration  procedure 
is  utilized.  The  spectral  peak  is  previously  known  and  provided  by  the  actual 
spectriim.  Thus, 

l=^^e'E(M  (B.7) 

To  start  the  nximerical  iteration,  initial  conditions  are  assigned,  by  experience,  as 

a  =  0.01,  0.01  <aa.  <  4,  0.01  <  at  <  4  (B.8) 

The  initial  value  -7  is  evaluated  by  Eq.(B.7).  Then,  a^  and  ai,  are  calculated  from 
Eqs.(B.5)  and  (B.6),  respectively,  utilizing  the  Bisection  method.  The  Bisection 
method  is  a  procedure  which  is  used  to  approximate  a  root  of  the  function  F{x]  in 
a  known  interval.  Suppose  that  F{x)  has  a  unique  root  between  Xi  and  X2,  there 
must  exist  the  relationship  F{xi)F{x2)  <  0.  Let  Xg  =  (ii+X2)/2;  if  F(xi)F(x<,)  <  0, 
there  must  be  a  root  between  Xi  and  x^.  Otherwise,  if  F{xi)F{xo)  >  0,  there  must 
be  a  root  between  xj  and  Xg.  Thus,  one  may  reiterate  the  procedure  imtil  the 
resulting  root  converges  in  a  satisfactorily  small  interval. 

By  using  Ca  and  a^  soved  from  the  Bisection  method,  two  solutions  of  a  can 
be  obtained  from  Eqs.(B.3)  and  (B.4),  respectively.  Taking  the  average  of  the  two 
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values  of  a  and  replacing  this  average  value  into  Eq.(B.7),  one  may  evaluate  a  new 
7.  Therefore,  one  may  reiterate  the  entire  procedure  until  an  acceptable  closure 
error  is  obtained. 


APPENDIX  C 
SUMMARY  OF  COMPUTED  AND  MEASURED  NON-DIRECTIONAL 

FREQUENCY  SPECTRA  AT  THE  BUOY  #41006  LOCATION 


C.l     Introduction 

The  wave  spectra  that  summarized  here  include  the  one-dimensional  frequency 
spectra  computed  by  the  UCWP  deepwater  sub-model  and  these  measured  from 
the  field  at  the  NOAA  buoy  #41006  location  in  Jajiuary,  September,  and  October 
1984.  Both  the  computed  and  measured  spectra  are  presented  by  graphs.  They 
are  managed  in  a  single  figure  if  occur  at  the  same  time.  The  eight-digit  number 
shown  in  esich  figure  indicates  the  time  when  spectra  occurred,  and  corresponds  to 
the  year,  month,  day,  and  hour,  each  represented  by  two  digits,  respectively.  The 
further  informations  shown  in  each  figure  are  the  significant  wave  heights  and  modal 
periods  corresponding  to  the  computed  and  measured  spectra,  respectively. 

C.2     Summary  of  Figures 

The  one-dimensional  frequency  spectra,  which  are  computed  and  measured  dur- 
ing the  period  of  January,  September,  and  October  1984,  are  summarized  in  the 
following  pages. 


128 


129 


NOBfl'UlOOe.     —    UCWP   MODEL 


LI  8110101 

13 

MS 

TN 

no- 

H.O 

n. 

1 

<4.2 

12. 

S 

rvi 

20- 

/; 

X 

'A 

V 

l\ 

a: 

1    \ 

X 

0- 

1      ^\^ 

sa^j 

i: 

»— 

L: 64010201 

(O 

HS       TN 

UJ 

HO- 

*t.i    9.1 

o 

H.%    J2.S 

>- 

o 

'» 

cc 

20- 

f  f 

LU 

*  t 

z 

^ 

LU 

/  \ 

0. 

/      -N.—-- 

L: 84010219 

UJ 

a. 

HO- 
20- 

HS       TN 
"S.t    11. I 
~S.«    I2.S 

0. 

/Vn^ 

I      '      1  ■  ' 

LlSUOlOSOll 

MS 

~2.9 

TN 
11.1 

2.7 

12.5 

K 

r 

eo- 

I.ia«0109l9 

Ht      TN 

10- 

~2.l     11.1 

~2.0    12.5 

n- 

n 

^N-r 

I.18II0I02IS 

MS      TN 
~3.S    11. 1 
~3.8    12.5 


Jk 


L:SU0102n( 

MS 
~3.3 
~3.S 

TM 

11.1 

12.5 

i^ 

L:8<I0]0219 

HS      TH 
~3.2    11. 1 
~3.3    12.S 


ik 


I.:e<t010221 

MS     TN 
~9.S    11.1 
"3.1    12.5 


Ltg<i010229 

HS       TN 
~3.1    11. 1 
~2.9    12. S 

'     /^ 

Ll  811010303 

NS       TN 

~3.1    It.l 
~2.5    12.5 

/^ 

LI8H01030S 

HS       TN 

~2.8    10.0 
~2.3    12.5 


Ll  811010307 

HS     TM 
~2.9    11.1 
~2.1    12.5 


Ll 84010308 1 

MS 

~2.8 

2.0 

TN 
11.1 

12. S 

-^^^ 

Li  84010311 

HS      TH 
~2.3    11.1 
~2.0    12.5 


-r 


Ll840103ts| 

MS 

TN 

2.8 

It.l 
12.S 

fK 

LlMOlOSnl 

HS 

TN 

2.8 
1.8 

12.5 

--^^ 

LlS40I03ia 

HS      TH 
~2.6    11.) 
~l.7    I2.S 


-^>^ 


Ll  84010321 

HS     TH 
~2.8    11.1 
"1.8    12.5 


0         0.1       0.2      0.3    0         0.1       0.2      0.3    0         0.1       0.2      0.3    0         0.1       0.2      0.3    0 

FREQUENCT  (HZ) 


j:v» 


LI  840 10323 

HS       TH 
~2.l|    11. I 
~i.7    12.5 


--^ 


0.1       0.2      0.3    0         0.1       0.2      0.3 


Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra. 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continuo). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 
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Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 


149 


Nafln«moo6. 


—    UCWP    MOOEL 


o 

<C      20 


20' 


Li  SO 1 02623 
HS     TH 

~2.2  a.i 
~i.»  11. 1 


Ll8«l027ll 

110- 

HJ      TN 
~.t    lO.O 
1.5    II. I 

20- 

U8«ia2723 

MS      TH 

110- 

1.5    10.0 
~2.«    1 1. 1 

X 

20- 

\ 

0. 

x: 

L:8«l02ail 

MS  TM 
~1.S  9.1 
~1.»    lO.O 


►— 

LI8II102823 

to- 

20- 

0. 

HS     TH 
~.3    8.3 
~.%    11. 1 

I 

u8iiia2aii 

HI  TH 
~.*  8.3 
~i.2    11. I 


U8«  102823 

HI      TO 

HO- 

~i.j  a.i 
~.o  n.i 

20- 

1        oil      o!2      0.' 

LI8XI027O1 

H9      TH 
~2.2    10.0 
1.7    11. 1 


JC^ 


■T- 


Ll8«102TI3 

H3      TM 
~1.7    10.0 
~2.5    11. 1 


-V>- 


T" 


Li SKI  02801 

MS      TH 

~.8    10.0 
~2.3    U.l 


-<-'^- 


T- 


L:8«I02813 

HS       TM 

~i.i»  9.1 

~1.7    10.0 


-r 


Ll  81(102901 

HS       TM 

~j.2   8.1 
~1.S    11.1 


Ll8«1029l3 
HS     TM 

~1.S  «.s 
~1.2    11.1 


T" 


Ul8«1030Ol 
H8     TH 

~i.»  a.i 

"1.0    11. I 


-r 


L!e<ii 

02703 

HS 
2.2 

"i.a 

TH 

10.0 

11.1 

h— 

L>mi027IS 
HS      TH 

~i.a  10.0 

~2.»    11.1 


LI S« 102803 

MS      TH 

~.S   9.1 

"2.3    11. I 

.:8m  02815 

HS      TH 
"i.S   9.1 

"i.a  10.0 


L:S<II02903 

HS      TH 

~1.2   9.1 

~1.3    11.1 

LI  80102015 

MS       TH 
""l.*    8.3 
"1.2    11.1 


LI8O103003 

HS  TH 
~1.3  7.1 
~i.O    11.1 


-r 


L:  am  02705 

HS      TM 
~2.1    10.0 
~2.8    11. I 


-'/y~ 


1- 


L184I027I7 

HS      TH 
~l.7    10.0 
~2.«    U.l 


i?i=. 


L18II10280S 

HS  TH 
~l.8  9.1 
~2.2    11.1 


.J4i^ 


-V 


L:8gi028n 

HS      TM 

~i.»  9.1 
~I.S   u.l 


-r 


L!aM10290S 

HS       TM 

~1.2   9.1 
~1.3    U.l 


"T" 


L180102S17 

HS  TM 
~1.3  9.1 
~i.l    U.l 


*r 


T- 


L1801030OS 

HS  TM 
~I.S  7.1 
~i.O  11.1 


L: 8*102707 

HS   TM 
~2.1  10.0 
~2.S  U.l 


Li  811102719 

HS   TM 
~i.7  10.0 
~2.»  U.l 


LI  80102807 

MS  TH 
~i.8  9.1 
~2.1  10.0 


■T- 


L: 80102819 

HS  TM 
~1.3  9.1 
~"l.S  U.l 


-r 


LI8O102907 

HS  TH 
~1.2  9.1 
~1.2  U.l 


T- 


Li 80102919 

MS   TH 

~1.3  9.3 
~l.l  U.l 


LI  80103007 

HS  TH 
~1.3  9.1 
~i.O  0.0 


T" 


0.2   0.3  0    0.1   0.2   0.3  0    0.1   0.2   0.3  0    0.1   0.2   0.3  0    0.1   0.2   0.3  0 

FREOUENCT  (HZ) 


L: 80102709 

HS   TH 
~2.0  10.0 
~2.S  U.l 


J4^ 


T" 


L: 80102721 

HS   TM 
~1.7  lO.O 
~2.0  U.l 


*T— 


LI  90102809 

HS   TH 

"l.S  9.1 
~2.0  10.0 


T- 


1:90102821 

HS  TH 
~1.0  9.1 
""l .  0  U.l 


-r 


Ll 90102909 

HS  TN 
"■|.3  9.1 
~i.2  U.l 


LI  80102921 

H9  TN 
~1.3  7.7 
~i .  1  U.l 


T" 


L 190103008 

HS  TH 
~1.«  9.9 
~i.O  0.0 


Figure  C.l:  Summary  of  computed  and  measured  frequency  spectra(continue). 


APPENDIX  D 
PLOTS  OF  WAVE  ENERGY  FLUX  ELEMENTS 


D.l     Summary  of  Figures 

The  rates  of  change  of  individual  energy  sources  on  the  frequency  domain  are 
shown  in  the  following  figures  for  the  six  csises,  Case  1  to  Case  6,  studied  in  Section 
4.5.  They  are  plotted  here  baised  on  the  result  of  numerical  hindccists  from  the 
UCWP  model.  Case  1  to  Case  3  are  in  the  deep  water  at  the  buoy  #41006  location. 
Case  4  to  Ceise  6  are  in  the  shallow  water  at  the  CDN  Marineland  station.  Case 
1  occurred  during  a  growing  period  under  increasing  wind  condition  at  1:00  a.m. 
Sep.  8,  1984,  when  winds  were  increasing  from  26  knots  to  31  knots.  Caise  2  is  a 
decaying  case  at  5:00  a.m.  Sep.  8,  1984,  while  onshore  winds  decreased  from  30 
knots  to  28  knots  in  two  hours.  Case  3  is  a  local  wind  wave  growth  case  at  23:00 
p.m.  Oct.  13,  1984,  when  winds  remained  a  constant  speed  at  12  knots  for  the  last 
twelve  hours.  Case  4  occurred  during  a  growing  period  under  increasing  winds  at 
5:00  a.m.  Sep.  8,  1984,  when  the  onshore  winds  increasing  from  27  knots  to  30 
knots  in  two  hours.  Case  5  is  a  decaying  case  at  1:00  a.m.  Sep.  10,  1984,  while 
onshore  winds  decreased  from  34  knots  to  29  knots  in  two  hours.  Case  6  occurred 
at  3:00  a.m.  Sep.  10,  1984,  when  onshore  winds  kept  dropping  from  29  knots  to  26 
knots  in  two  hours. 
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Figure  D.l:  Plot  of  deepwater  energy  fluxes  in  growing  period;  Ccise  1. 
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Figure  D.2:  Plot  of  deepwater  energy  fluxes  in  decaying  period;  Case  2. 
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Figure  D.3:  Plot  of  deepwater  energy  fluxes  during  small  winds;  Case  3. 
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Figure  D.4:   Plot  of  individual  shallow  water  effect  sources  during  growing  period; 
Case  4. 
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Figure  D.5:  Plot  of  individual  shallow  water  effect  sources  during  decaying  period; 
Case  5. 
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Figure  D.6:  Plot  of  individual  shallow  water  effect  sources  during  decaying  period; 
Case  6. 
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